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CROSS SECTION FOR THE REACTION Al?(y, n)AP*? 


By R. N. H. HaAstam, W. N. RoBerts,? anp D. S. Ross? 


ABSTRACT 


The cross section for the reaction Al*"(y, 2)Al** has been remeasured by the 
determination of Al** activity. The results are in agreement with an earlier value 
obtained by the activity method, but differ from the result obtained by neutron 
detection. The discrepancy is thought to be due to complexity in the Al?* decay 
scheme. The reaction threshold is measured as 13.4 + 0.2 Mev., the maximum 
energy of positrons from Al* is 3.2 + 0.1 Mev., and the half-life is determined 


as 6.5 + 0.1 sec. 
INTRODUCTION 


Cross sections for gamma-neutron reactions may be measured either by 
detection of the emitted neutrons or by determination of the residual activity. 
In the several cases in which both methods may be applied, satisfactory 
agreement is generally obtained between the results of the two methods (11). 
An exception is the reaction Al?’(y, ”)Al**, for which the neutron detection 


method gives values of the peak and integrated cross sections nearly three 
times those obtained by the residual activity method (7, 11). This discrep- 
ancy is far greater than the estimated probable errors of the two experiments 
should allow. In view of the theoretical implications of this difference, 
should it prove real (13), it was thought worth while to repeat the deter- 
mination of the cross section of Al®’(y, )Al** by the residual activity method. 


EXPERIMENTAL AND RESULTS 


A relative activation curve for the reaction Al®’(y, ”)Al** was obtained by 
irradiation of the pure aluminum samples at maximum betatron energies 
ranging from 13 to 25 Mev., and counting of the resulting short-lived positron 
activity corresponding to each energy setting. Owing to the relatively short 
half-life of the Al* (6.5 sec.), special apparatus is required to determine the 
activity. 

An aluminum disk, mounted at the end of a light swinging arm, actuated 
by an electromagnet, is alternately raised into the betatron beam and lowered 
into position above a carefully shielded thin-window counter. In a typical 
cycle of events the disk is irradiated from 0-15 sec., lowered into position, 
and counted from 18-30 sec. by one scaler and again from 33-57 sec. by a 
second scaler. The second scaler follows any variation in background due to 

1Manuscript received February 18, 1984. 

Contribution from the Physics Department, University of Saskatchewan, Saskatoon, Sask. 
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longer period activities. Both scalers use the same preamplifier circuit. The 
disk is raised into position for a second irradiation by the end of 60 sec. This 
whole cycle is repeated as many times as necessary to secure good counting 
statistics, and is controlled by a switching mechanism driven by a synchronous 
motor. 

During irradiation the aluminum disk is shielded with paraffin wax and 
cadmium to reduce the 2.3 min. activity due to Al”"(, ~)Al*. The thickness 
of the disk (450 mgm./cm.?) is chosen to give a good counting rate without 
too large a background due to annihilation radiation. The X-ray dose given to 
the sample is measured by a Baldwin-—Farmer ionization meter.* The betatron 
injection voltage, and thus X-ray yield, is on only during the irradiation period. 
At a fixed betatron operating energy the accumulated counts registered by 
scaler 1, corrected for background, and divided by X-ray dose, is remarkably 
constant. 

Relative activation curves were obtained also by the use of two other types 
of apparatus. These are modifications of apparatus previously described 
(4, 5). One of these is a “‘dropping’’ apparatus, in which an aluminum cylinder 
is alternately raised into the betatron beam and dropped into counting position 
over a thin-walled cylindrical Geiger tube. The cycle may be repeated as 
often as desired to secure good counting statistics. This apparatus proved 
particularly suitable for threshold determination owing to the good counting 
geometry, but did not lend itself to the determination of the positron energy 
or the absolute disintegration rate of Al**. The shapes of the relative activation 
curves obtained with all three types of apparatus were identical, within the 
limits of experimental error. 

The relative activation curve was normalized to the activity resulting from 
Cu®(y, 2)Cu®?, which is used as a standard in this laboratory. This was done 
by irradiation of disks of copper and of aluminum of the same diameter at 
the same betatron energy (20 Mev.) and at the same position. The resulting 
positron activities were counted in the same geometry. Corrections for self- 
absorption were applied using the method of Baker and Katz (2). This 
eliminated the necessity for the determination of a self-absorption curve in 
copper and in aluminum, as was done by Katz and Cameron (7) in their 
study of this reaction. Self-scattering corrections cancel out if the sample 
thickness is greater than 0.2 times the positron range in each case. 

In the application of the normalization procedure it is necessary to know 
the half-life of Al** and the maximum positron energy. The half-life was 
measured as 6.5 + 0.1 sec., in excellent agreement with the value given 
by Katz and Cameron (7), but somewhat lower than the recently reported 
value of 6.68 + 0.11 sec. (3). 

Previously reported values of the positron energy have covered a wide 
range (1). In the present experiment the energy was measured by determina- 
tion of an aluminum transmission curve with the “arm’’ apparatus, filters 
of various thickness being interposed between the irradiated disk and the 


*Product of the Baldwin Instrument Company, Dartford, Kent, England. 





HASLAM ET AL.: Al*"(7, n)Al’® CROSS SECTION 363 


counter window. Analysis of the transmission curve by the method of Katz 
and Penfold (8) yielded an end-point energy of 3.2 + 0.1 Mev. for the Al* 
positron. 

The reaction threshold was measured by use of the dropping apparatus, 
as 13.4 + 0.2 Mev. This value is considerably higher than the value 12.75 
+ 0.20 Mev. obtained by Sher et al. (12) using neutron detection. McElhinney 
et al. (10) obtained a threshold of 14.0 + 0.4 Mev. from observation of the 
Al® activity. Threshold measurement by the activity method is complicated 
by a considerable background, due probably to Al*’(m, y)Al*. However, a 
threshold of 13.4 + 0.2 and a positron energy of 3.2 + 0.1 Mev. are con- 
sistent with the masses given by Li (9) for Al?” and Mg”, which predict a 
difference of 10.1 Mev. between threshold energy and positron kinetic energy. 

Analysis of the activation curve by the photon difference method (6) leads 
to a cross section of the usual resonancelike shape, whose characteristics 
are compared in the following table with the corresponding values given by 
other investigators. 


Maximum Position Integrated - 
cross of cross alf- 
Method Reference section maximum section width 
(mbarns) (Mev.) (Mev-barns) (Mev.) 
Residual activity Present work 8.4 19.5 0.042 5 
Residual activity (7) 8.1 19.2 0.045 4.7 


Neutron detection (11) | 24.4 19.7 0.132 4 





DISCUSSION 


It is seen that the results of the present experiment are in substantial 
agreement with those obtained, also by the residual activity method, by 
Katz and Cameron. Maximum and integrated cross sections are about one- 
third as large as the corresponding quantities obtained by neutron detection. 

The integrated cross section, measured by neutron detection, may contain 
a small contribution from the (y, mp) reaction. Since this reaction has a 
threshold of about 19.5 Mev., it could not appreciably affect the maximum 
cross section value. 

The details of the present experiment were checked to seek a possible 
explanation for the smaller values. At all times counting rates were kept so 
low that counting losses were negligible. The aluminum and copper activities 
were counted on the same day with the same equipment. It is true that the 
values of the Al” half-life and maximum positron energy enter into the 
calculation of absolute activity. A decrease of 1 Mev. in assumed positron 

-energy would lead to a 40% increase in activation value, and thus in cross- 
section value; a decrease of 0.5 sec. in the assumed half-life would lead to a 
7% increase in cross section. These are thought to be extreme possible varia- 
tions, and cannot account for the difference between the results obtained by 
the two methods of measurement. 
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Assuming the validity of the neutron detection results, which have been 
very carefully checked, it appears that the discrepancy must be explained 
by the assumption that the decay scheme of Al* is not simple. This is in 
agreement with the suggestions recentiy advanced by Stahelin (13). 


The authors would like to thank Dr. L. Katz for valuable discussions. They 
are indebted to the National Research Council of Canada for financial support. 
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ON THE THEORY OF AN ANTENNA WITH AN INFINITE CORNER 
REFLECTOR! 


By JAMEs R. Wait 


ABSTRACT 


Some comments are made concerning the theory of a dipole radiator inside 
a corner or ‘‘V”’ reflector of infinite extent. A straightforward solution is obtained 
for the resultant fields anywhere within the angle subtended by the V. The 
result is discussed and compared with Moullin’s theoretical treatment for the 
filament source. 


INTRODUCTION 


To obtain a simple directive radiating system for operation in the high- 
and very-high-frequency radio ranges, a corner reflector antenna is sometimes 
employed. The antenna is mounted on the bisector of two flat sheets meeting 
in a “V’’. Kraus (1) and Moullin (2) discuss this type of antenna and employ 
the theory of images to calculate field patterns that would actually only 
apply to the case where the flat sheets were of infinite extent. These ideal 
solutions form a valuable guide to the behavior of the corner reflector antenna 
of finite size. As pointed out by Moullin, the image method appears to be 
applicable only for angles that are a submultiple of x. He is rather skeptical 
as to whether the results expressed as a series of Bessel functions can be 
generalized to corner reflector angles of an arbitrary value between 0 and 
27. It is shown in this note that the theory can be extended to any angle by 
simply admitting multivalued solutions of the wave equations in cylindrical 
coordinates, at the outset. In addition Moullin’s restriction to a line source 
excitation is removed and the antenna is represented instead by an electric 
dipole oriented in a direction parallel to the apex of the reflector. The resulting 
formal solution which contains an infinite series of integrals is evaluated for 
the far field by applying the saddle point method of integration. The final 
result is then expressed as a singly infinite series containing a Bessel function 
which in general has a fractional order. For antennas that are within a few 
wave lengths of the apex the series converges very rapidly. Several repre- 
sentative pattern calculations are illustrated. 

A similar problem has been treated by Tilston (3) who employs spherical 
wave functions and expresses his final result as a doubly infinite series con- 
taining fractional order Ferrer’s Associated Legendre functions and Bessel 
functions. He gave no numerical results. 


THE FORMAL SOLUTION 


With respect to a cylindrical coordinate system (p, ¢, 2) the corner reflector 
is defined by the planes ¢ = 0 and ¢ = y which intersect along the z axis. 
The dipole of electric moment # is situated at the point C(po, $0, 0), where 
¥ > do > 0, and is oriented in the z direction. The primary Hertz vector of 
the dipole at the point P(p, ¢, z) has only a z component II, and is given by 


1Manuscript received February 3, 1954. 
Contribution from the Radio Physics Laboratory of the Defence Research Board, Ottawa, Canada. 
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[1] Il, = pe-*"/r 
where r = (p:2 + 22)4 
and p1 = [p? + po? — 2pop cos (go — ¢)]?. 


To satisfy the boundary conditions it is only necessary to assume that the 
secondary field can be represented by a Hertz vector which has also a sole 
component in the z direction and is denoted by II;. The resultant Hertz 
vector, II = II, + Il,, is a solution of the homogeneous wave equation except 
at the dipole and therefore, in general, is made up of combinations of 
+igz COS 

[2] Z,[uple"™" -. 6 
where u = (8? — g?)* and Z, is a cylindrical Bessel function of order v. 
Since II is not necessarily periodic in ¢, v is not restricted to an integer. For 
a perfectly conducting flat sheet the boundary condition for this scalar prob- 
lem is simply 

Il = Oat ¢ = Oand y. 


The angle ¢ then enters the solution as sin v@ where vy = mz/y and m is any 
integer. The orders of the Bessel functions are thus specified. The primary 
field is now written as an integral representation 


3) — 


where H®) is the Hankel function of the second type of order zero. The in- 
tegration is along the real axes of g and @ and is assumed to be slightly com- 
plex to avoid difficulty with branch points. It is now convenient to employ 
an addition theorem for the Hankel function given by 


$e —igz 
i Ho (up) e dg 


—™<co 





f 


[4] H@ (up:) = 7 x im H(upo) J,(up) cos v(¢ — ¢o) 
for p < po where €9 = 1, €m = 2 for m = 1, 2, 3,....The arguments upp and 


up are interchanged if p > po. This result is well known (5) if v is an integer 
(i.e. Y = m) and it can be readily checked for arbitrary v by noting that the 
right-hand side: 
1. is a solution of the wave equation, 
2. is continuous at p = po for d ¥ do, 
3. is symmetrical in ¢ about @o, and 
4. has the proper singularity for p — po and ¢ > @o. 
The primary field is then written 
. oo +00 
i a 
[5] = — FEY eu cosv(o — bs) J HP (ups) So(up) "dg 
for p < po and z > 0. The factor e+”? is employed if z < 0. 
The secondary field II, is now constructed so that ¢ appears as sin vd 


in the expression for II. It is evident that 
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[6] II, = ore. Em COS v(h + go) 7 H® (upo) J.(up) e~ ™ dg 


for p < po and z > 0. 
The resultant field for the region p > po is then given by 


+o 
[7] II = - Be z 4 sin v¢ sin vo f J,(upo) H? (up) e~ ** dg 


where u = (6? — and v = mx/y. The complete formal solution for the 
electric and magnetic fields, E and H, is then specified by carrying out the 
operations 
= (8? + grad div) II i,, 
H = (6?/nw) curl II i, 


where i, is a unit vector in the z direction. 


THE RADIATION FIELDS 


Since in the practical situation the fields are always to be observed at large 
distances from the antenna, it is possible to simplify the formal situation to 
a useful form by an approximate evaluation of the integrals. The integral 
to be evaluated is 


: +a 
[8] =—2] J,(upo) H® (up) e~*dg 


and since 8p > 1 it is permissible to use the approximation that 


: - oo ss 
H®? (up) = ai 2 gi rer. 
Tup 


The integral is then in a form to be evaluated by the saddle point method 
if the integration is transformed to the complex a plane defined by 
= B cosa. 


The integral J is then given by 
ive /2 1434/4 G1 - 28 sin a . —i8Reos(—a) 
I=e"e 2 —————- J, (Ap cin a)¢ da 
7—t2 Tp 


with p = R sin 6, = Ros 6, and where the contour is along a line parallel 
to the negative imaginary axis, from z to 0 on the real axis, and then along the 
positive imaginary axis. The integral has a saddle point at a = @ and the 
contour is now transformed to the path of steepest descent defined by 

cos (@— a) = 1 — ix? 


‘where x ranges from —© to +o through real values. It can be verified 
that no poles are crossed in this deformation of the contour so that 


if (ast sina)! J, (Bp sin poh, ARE dee. 


#6 ee Q—i/2)) 


The term in square brackets in the integrand can be expanded in a power series 
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in x? and the integration is then carried out term by term. Retaining only the 
leading term, the integral is given by 
[9] I =e" ¢ ** R™ J,(Bpsin 0) 


where the remaining terms contain higher powers of R-. 
Employing this result the Hertz vector at large distances from the antenna is 


[10] T= pre RG 
where 
[1 1] = : » sin — sin ws a /wy Iz (Bpo sin 9) 


and R = (p? + 2?)}, 6 = tan-'(p/z). The only field components that vary as 
1/R are given by 
8° all B’pr e #8 








[12] Ay = — ea ae = a ae sin@ G 
and 
[13] Eo Bpre**® R™' sin G. 


The shape of the radiation pattern is therefore characterized by the function 
G, which is readily computed since the series is rapidly convergent if po sin 6 
is within a few wave lengths. An interesting special case is when y = 27 
which corresponds to a knife edge where Bessel functions are of order m/2. 
This particular result could also be obtained by applying the reciprocity 
theorem to Sommerfeld’s famous problem of a plane wave incident on a 
knife edge. 

The function G has actually been given by Moullin who derived it by an 
application of the image theory. He considered cases where ma/y was an 
integer which corresponded to corner reflector angles of 15°, 30°, 45°, 60°, 
90°, and 180°. Also he considered this antenna to be represented by a filament 
source parallel to the apex of the reflector. As seen from equation [11] this 
case corresponds to G when @ = 90°. It is interesting to note that Moullin’s 
method of images can easily be extended to arbitrary angles of y if the con- 
cept of multifold space is admitted. For example in the case of a knife edge 
if the source dipole is at (po, 0, 0) then @o is taken in the first branch (0 <6) < 27) 
of a twofold Riemann space. The image dipole at (po, @o’, 0) is then taken 
in the second branch ( —2m < ¢o’ < 0) where ¢o’ = —¢o to satisfy the boundary 
condition on the screen. 

PATTERN CALCULATIONS 


The function G has been calculated for values of y (30°, 45°, 60°) which 
would correspond to typical corner reflector angles. Five different cases are 
considered: Bpo sin @ = 6, 12, 18, 24, and 29, which are henceforth designated 
by I, II, III, IV, and V respectively. The factor e***/™ is either + 7 or + 1 
in each of these series to be summed, which simplified the calculations con- 
siderably. The terms sin (mr@/y) sin (mago/W) depend only on ¢/y and hence 
the same values can be used for each y. Also since the dipole is taken to be 
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PATTERN FUNCTION G 









“O 5 10 15 
ANGLE 


Fic. 1. The pattern function G as a function of the azimuthal angle ¢ for a dipole located 
on the bisector of a corner reflector and parallel to the apex for a corner angle y of 30°. 
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Fic. 2. The pattern function G for a corner angle y of 45°. 
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PATTERN FUNCTION Gi 


Fic. 3. The pattern function G/i for a corner angle y of 60°. 


on the bisector (¢o = ¥/2), the function G is symmetrical about ¢9 so that 
the calculations are only carried out for values ¢ between 0 and @o. 

The function G is plotted in Figs. 1 and 2 as a function of ¢ for Y = 30° and 
45° respectively. The function G/z is plotted in Fig. 3 as a function of @ for 
y = 60°. In displaying the pattern in this form the information regarding 
the phase of the radiated field is not lost as would be the case if a conventional 
polar diagram were made. It is interesting to note that for the case I 
(i.e. Bp sin 8 = 6) the pattern is practically a simple sinusoid since only the 
first term in the series is significant. In fact this simple pattern is a satisfactory 
representation when po is less than about a wave length for reflector angles 
y less than 60°. 

The general behavior of the patterns is not surprising. As the distance 
between the dipole and apex is increased the pattern splits into lobes. Usually 
these secondary maxima are undesirable from the standpoint of conventional 
operation of directive antennas. However, they could possibly have some use 
in special applications such as direction finders where the center lobe can be 
made sharper by simply increasing the apex—dipole distance po in regular steps. 

Calculations for other angles of y could be readily carried out with the aid 
of a table of Bessel functions of both integral and fractional order. Examples 
of such computations have been given previously for a dipole near a con- 
ducting half-plane or knife edge (4). For small values of po (less than 4/4) 
the pattern is again given approximately by a sinusoid, G being simply pro- 
portional to (sin (r¢/yW) sin (4¢0/yP)). 
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The effect of employing finite sheets rather than infinite sheets has been 
discussed by Moullin (2), who carried out a series of interesting experiments 
for corner angles of 45° and 60°. Apparently the lobes beyond the edges of 
the sheet are only a few per cent of the central maximum if the sheets are 
greater than 5\/4 in height and 3\/2 in width. These spurious lobes become 
more significant for smaller corner angles and cases where the dipole-apex 
distance becomes greater than a wave length. 
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DIFFRACTION OF 3.2 CM. ELECTROMAGNETIC WAVES BY 
CYLINDRICAL OBJECTS! 


By S. T. WILEs? AND A. B. McLay 


ABSTRACT 


Diffraction patterns of a brass tube and a hard rubber rod, each a cylinder of 
1 in. diameter, in a nearly plane beam of square-wave modulated 3 cm. waves 
with electric vector parallel to the cylinder axis, have been measured in several 
planes transverse to the incident beam direction. Experimental results for the 
conducting cylinder agree closely with calculations based on scalar diffraction 
theory. Patterns of the dielectric rod show a pronounced central peak immedi- 
ately behind the rod and other intensity effects differing from the conducting 
cylinder patterns, particularly in the vicinity of the shadow. 


1. INTRODUCTION 


Development of ultra-high-frequency sources of electromagnetic waves with 
lengths in the centimeter range has made possible investigations paralleling 
those of optics but under experimental conditions not practicable when using 
light waves. In particular, diffraction effects are very pronounced and can be 
studied in the close vicinity of an object with one or more dimensions of the 
order of the incident wave length to within a distance of less than a wave 
length from its surface.* 

In the years since 1945 increasing attention has been paid to experimental 
and theoretical studies of microwave diffraction. Much of this has been 
concerned with diffraction by plane apertures and obstacles, particularly those 
with circular shape, and by an infinite half-plane. Experimental results have 
directed attention to known or suspected inadequacies of the Huygens- 
Kirchhoff theory of diffraction and to the necessity of a more exact treatment. 
Less attention has been paid to diffraction by three-dimensional objects and 
much of this has been confined to far-field effects because of their practical 
importance in radar wave propagation. Diffraction by long right-circular 
cylinders in incident plane radiation is the simplest to treat theoretically but 
computation of exact very-near fields will in general be extremely laborious, 
since approximation methods applicable to far-field effects are of little assis- 
tance. For the special case of a long cylinder in a field incident normal to and 
plane-polarized parallel to the long axis of the cylinder, two-dimensional scalar 
theory can be used and computation is considerably reduced. 

A program of investigation of diffraction of 3.2 cm. microwaves by solid 
objects, beginning with long right-circular cylinders with diameters of the 
order of the wave length, is being undertaken in this laboratory. In this paper 
results are reported on the diffraction fields of a brass conducting cylinder and 

1 Manuscript received January 25, 1954. 
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*Such very-near-field diffraction can also be studied using acoustic waves with lengths comparable 
to those of electromagnetic microwaves, but the vector nature of the laiter generally makes inter preta- 
tion of the observed effects much more difficult. 
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of a dielectric cylinder of hard rubber, both 1 in. in diameter. Calculation of 
the former field has been made and compared with the observed field. Know- 
ledge of its theoretical field affords a very valuable check on errors introduced 
by the many disturbing effects met with in experimental studies of microwave 
diffraction. 

Results of an investigation closely paralleling the present one have been 
reported by Kodis (3) who also used the calculated field of a conducting cylin- 
der to test the performance of his diffraction apparatus, which was set up on 
the half-plane principle for use with 1.25 cm. radiation. He also made some 
observations with polystyrene cylinders and with a conducting infinite half- 
plane. His experimental arrangement is a better one than ours for the particular 
case when the incident radiation is polarized parallel to the cylinder axis. It is 
not so easily adaptable to study of 3.2 cm. waves. With the arrangement de- 
scribed below any polarization of the incident radiation may be used and also 
the field can be observed in the plane transverse to the incident radiation that 
contains the cylinder axis, as well as in planes behind the cylinder. Knowledge 
of the simple patterns in the transverse in-axis plane of a dielectric cylinder may 
be useful for interpretation of the more complex patterns in other planes and 
computation of some values in the simple pattern could be a useful check on 
experimental accuracy. 


2. APPARATUS AND EXPERIMENTAL PROCEDURE 


A 723A/B reflex klystron oscillator, with a TVN-7BL power supply and 
square-wave modulator, built from specifications by Montgomery (4) was 
used to generate 9375 Mc./s. oscillations. These were modulated at 1000 c./s. 
The 3.2 cm. waves were radiated by a rectangular H-plane brass horn, 80 cm. 
long from mouth to geometrical apex, flared at 20°, and mounted with the 
H-plane horizontal. Thus the axis of propagation was horizontal and the field 
vertically plane polarized. A cavity absorption wavemeter was used to set the 
frequency. A number of different tests showed the oscillator stability to be 
better than 1 part in 10,000 over a much wider range of conditions than those 
in the actual experiments. 

A IN28A crystal diode was used as a probe to measure relative field inten- 
sities directly, considering it to be a square-law detector. The brass at one 
end of the crystal cartridge was turned down to form a rod to match the rod 
terminal at the other end, as an approach to a simple dipole antenna. The two 
end rods were then cut by trial to obtain best experimental matching to coaxial 
cable. The crystal was mounted in a slot in the wedge-shaped end of a hard 
rubber tube 75 cm. long, and the cable was led from the probe through the tube 
away from the source of radiation, then down to the floor and to the input of a 
broad band a-c. amplifier. The amplifier has maximum gain of 35,000. Its 
output was clamped by a diode with variable d-c. bias, adjusted to balance 
out the rectified noise of the amplifier at given gain setting. A small part of 
the output was used to operate a Brown Electronik Potentiometer, Type 
153X12V-X-30 with 5 mv. range. 

The probe assembly was mounted on a 2 in. X 2 in. wooden stand, coated 
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with aquadag, and set in a heavy optical bench base. A 3000 r.p.m. induction 
motor, a 9900 : 1 reduction gear, a hard rubber drum, and a stranded string 
served to drive the probe parallel to the bench track at about 1 mm. per sec. 
The string was wound a few times around the drum and over a pulley at the 
other end of the track. Tension was obtained and held constant by use of a 
spring balance attached to the pulley. The relative speed of chart and probe 
was nearly but not quite constant throughout many experimental runs. Cali- 
bration marks were therefore put on each chart record at a synchronized 
starting time of probe and chart and at a later time of run. The track was 
placed on the floor so that the probe moved in a plane at right angles to the 
axis of propagation of the horn radiator and at 615 cm. (about 192d) from its 
geometrical apex. This plane of motion was kept fixed for all diffraction runs 
and the distance from the long axis of the vertical cylindrical diffracting object 
to the plane was varied by moving the object. The minimum distance from the 
surface of an object to the center of the probe was limited by the radius of the 
probe cartridge to about 4 mm. 

Before setting up for final intensity runs, the phase in the radiation pattern 
of the horn was investigated by a simple method described by Montgomery 
(5). A simple half-wave dipole of copper wire was matched to coaxial cable by 
a choke coupling of slotted thin cylindrical copper tube. The dipole length and 
corresponding matched slot length were obtained by cut and try experiments. 
A study of the phase fronts at several distances from the horn showed that these 
were sensibly cylindrical with center at or near the apex of the horn, at least 
over a wider angle of the beam than that used in making intensity measure- 
ments. Kodis (3) obtained the same result by more exact phase measurements 
at a distance 80 \ from a pyramidal horn. He also measured phases as well as 
amplitudes in the diffraction patterns and found close agreement with cal- 
culated values for conducting cylinders if the source was considered to be at the 
apex of the horn. We did not have equipment suitable for accurately measuring 
the rapidly changing phase in our diffraction patterns but comparison of ob- 
served and calculated intensities supports the same conclusion as to source 
position. 

Preliminary trials of the apparatus were made in a laboratory measuring 
20 ft. by 20 ft. but the equipment was moved to a larger room 25 ft. wide and 
75 ft. long for the final experiments. Effects of scattering and reflection from 
the walls and from other fixed objects near the walls and ceiling of this room 
were found to cause negligible disturbance. 

Intensity distributions in the diffraction fields of a cylindrical brass tube with 
1/16 in. walls, taken as the equivalent of a solid rod, and of a solid cylindrical 
rod of hard rubber, both 6 ft. long and of 1 in. outside diameter, have been 
obtained, with incident radiation polarized parallel to the long axis of the 
cylinders. Each cylinder was placed with its axis vertically intersecting the 
axis of propagation of the horn. Runs were made with the cylinder axis in the 
plane of motion of the probe and at selected points in front of this plane, i.e., 
towards the source. The positions were more closely spaced in the case of the 
dielectric rod because of the rapidly changing pattern in its vicinity. 
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After each diffraction run the cylinder was moved away and both probe and 
chart moved back to the original starting point. A run was then made to 
obtain the intensity of the incident field. This was quite constant for all runs 
and fairly uniform over more than 50 cm. but fell off outside this range closer 
to one side of center than to the other. Also the intensity showed small irregular 
variations which could not be eliminated. The runs were made over the range 
for which the incident field was most uniform from 15 cm. on the left to some 
40 cm. on the right. Runs in the plane of the cylinder axis were made on the 
extended right side only. Both sides should have been investigated in this 
case because of a noticeable asymmetry found in all the other patterns. 

Each of the runs with and without the rod was repeated at least once. All 
the runs were performed one after the other as quickly as possible to keep 
conditions closely the same in all. Immediately after completion of these the 
apparatus had to be moved from the large room. A rerun of the whole experi- 
ment with certain altered conditions would have been advisable but was not 
possible in that fairly ideal location. 

The dielectric constant, K,, and loss tangent, tan 6, of a small piece of the 
hard rubber rod were measured. The values found were K, = 3.45 and 
tan 6 = 0.104. The method of Dakin and Works (2) was used, which is 
claimed to be satisfactory for materials with losses less than or up to about a 
value tan 6 = 0.1. 


3. CALCULATIONS AND EXPERIMENTAL RESULTS 


The geometry of a diffraction experiment in the horizontal H-plane of the 
horn radiator is shown, not to scale, in Fig. 1. S is the source of 3.2 cm. radi- 





Fic. 1. Geometry of a diffraction experiment in the H-plane of the horn radiator. 


ation, taken to be at the geometrical apex of the horn. SX is the axis of propa- 
gation outward along the projected geometrical axis of symmetry of the 
horn. The probe P was moved at constant speed transverse to SX at distance 
Xo from S, which was fixed for all runs. The diffracting cylinder was placed 
with axis O on SX at selected small distances x from the fixed line of motion of 
' P. Ois taken as origin of rectangular and cylindrical coordinates with common 
Z axis normal to the plane of the diagram. It is more convenient experimen- 
tally and for purposes of calculation to vary x by moving the cylinder than by 
moving the track and probe assembly with the cylinder in a fixed position. 
The waves from S are considered to be cylindrical at P, shown by the broken 
arc with sagitta s, for the range of positions of P actually used, but plane to a 
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close approximation over the cross section of the cylinder because of the small 
value of a, 1.27 cm., and the small values of x used in the experiment. The 
incident waves are linearly polarized parallel to the long axis of the cylinder, 
the Z axis. The incident amplitude (modulus) is taken to be constant at all 
positions of the probe P since s < Xo and at all points on the surface of the 
cylinder because x < Xo and a is small. The waves scattered by the cylinder 
will also be linearly polarized parallel to the Z axis. The scattering problem 
thus reduces to a two-dimensional one and is solvable by scalar field theory. 
The theory of scattering of plane electromagnetic waves by an infinitely 
long, perfectly conducting circular cylinder with axis parallel to the electric 
vector of the incident field can be obtained using general principles given by 
Stratton (7). Or a solution given by Morse (6) for acoustic scattering may be 
altered to suit the electromagnetic case by suitable substitution of electro- 
magnetic magnitudes (vector components) and boundary conditions for 
acoustic ones. The space dependent part of a plane incident wave at 
P(x, y;7r,0) with modulus Ey, taken as unity can be replaced by a set of 
cylindrical waves using cylindrical coordinates about the OZ axis such that 


. ikz ikr cos @ 
Eine = € <0 


> i” In (kr) cos mb 


m=—co 


(1] 


where J,,(kr) are Bessel functions and k = 27/X. For the wave at P scattered 
by the cylinder one takes the general set of outgoing waves given by 


[2] Ese = > Am Hm (kr) cos mb 
m=0 


where A, are amplitudes and H,,“ (kr) is the Hankel function of the first 
kind. The resultant field E at P is: 


[3] E= Eine + Eve: 


Since Ey, and E,, are both tangential at all points on the surface of the 
cylinder, 


[4] E=0Oatr=a. 

By solution of [1], [2], [3], and [4] 

. — — Jolka)_ — — 9;n _Jm(ka) 
[5] Ao = H.” (ka) and Ay= 21 H,,” (ka) ; 


If now the incident wave is taken to be cylindrical at P rather than plane as 
above, 
[6] Eine = e*?+® where s = y?/2Xo if Xo is large, 


but if it is sufficiently plane over the cross section of the cylinder, the scattered 
wave can still be represented by [2]. Thus: 


[7] E= eto _ SS ™ z — H,," (kr) cos mé. 
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By use of this equation the value of the modulus | E] and phase of E at P can be 
determined. 

Computation of fields using Equation [7] is very laborious, but necessary for 
evaluation very close to the cylinder since approximate solutions will not 
suffice when kr is small. Convergence was sufficiently complete using order 
numbers of functions up to m = 7 for arguments ka = 2.494 and kr < 5 and 
to m = 6 for larger kr values, to keep well below experimental error. Values of 
the functions for arguments up to 25 were obtained from the British Associa- 
tion Mathematical Tables (1). 

Calculated intensity distributions in the diffraction patterns of a conducting 
cylinder with radius a = 1.27 cm., ka = 2.494, are shown to values kr = 25 
by the curves in Fig. 2 for three values of x, namely, x = 0, 1.67 cm., and 
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Fic. 2. Comparison of experimental and calculated intensities in diffraction patterns of a 
brass cylinder. 


7.67 cm. Values of |E 2 calculated from Equation [7] are the ordinates, but are 
expressed as the ratios J/J, of the intensity at a point in the diffraction field 
to the intensity in the incident field at the same point to stress experimental 
measurements in which J; was not exactly constant. 

The experimental measurements of patterns of the conducting cylinder, 
corresponding to the calculated ones, are shown in Fig. 2 by crosses and dots 
for the two sides L (left) and R (right) of center respectively. Experimentally 
I and J, are ordinates of a given curve at each displacement as measured 
directly from the recorder chart. These and all other experimental values are 
averages obtained from measurements of two chart records in each case. The 
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range of ky values in Fig. 2 is sufficient to include all measurements made on 
the left side. Measurements made out to much larger ky values on the right 
sides showed that peaks (and troughs) in any pattern were of nearly uniform 
intensity and nearly equally spaced. Because of the observed asymmetry closer 
to the cylinder, the spacing is probably a little larger than predicted. Records 
were also made at two larger distances x — a equal to 5X and 10X. These 
show the expected increase of spacing of maxima and small increasing inten- 
sity in the immediate shadow of the cylinder. Since results for the conducting 
cylinder are not claimed as original but are given only as a test of performance 
of the equipment, measurements for values of kr larger than 25 have not been 
reported here. 

Experimental results of diffraction measurements in the case of the hard 
rubber cylinder (corresponding to those in Fig. 2) are shown in Fig. 3. The 








Fic. 3. Experimental intensities in diffraction patterns of acylinder of hard rubber with 
K, = 3.45, tan é = 0.104. 


same kind of right—left asymmetry is evident. To calculate the patterns would 
require introduction of a third wave, a transmitted one, in addition to the 
incident and scattered waves, with account being taken of the relative di- 
electric constant, K, = 3.45, and loss tangent, tan 6 = 0.104. Computation of 
the resultant field from an equation resulting from solution of the three-wave 
case would be exceedingly tedious and time-consuming. Experimental results 
in such a complex case can be obtained much more quickly and easily than 
theoretical ones. In addition to the results in Fig. 3, measurements on the 
right sides only of all observed patterns out to y = 30cm. (ky = 60) are shown 
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in Fig. 4 to illustrate interesting changes in pattern with increasing x. Some 
of these details would not appear if the patterns were plotted to ky = 25 only, 
as in Fig. 3. The patterns are presumably all spread out a little because of 
asymmetry, which shows in the close-in parts of three of them in Fig. 3. 
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Fic. 4._ Experimental intensities in the right side only of diffraction patterns of a hard 
rubber cylinder. 


4. DISCUSSION OF RESULTS 


Inspection of Figs. 2 and 3 shows clearly the asymmetry of all these experi- 
mental curves, which no doubt resulted from an asymmetry of the field of the 
horn about its geometrical axis of propagation. If the left-right results for the 
conducting cylinder were averaged, the resulting curves would have peak and 
trough displacements, ky, close to the corresponding ones of the calculated 
theoretical patterns shown in Fig. 2. However, the experimental intensities of 
.all maxima are a little higher and those of minima lower than the calculated 
ones. Deviations of the same order are also noticeable in closely comparable 
results of Kodis (3). These errors could be caused by a small deviation from 
square-law response of the crystal rectifier to the very small input. The 
exceptional height of the first maximum in the case x = 0 probably indicates 
reradiation between probe and surface of the cylinder, which are closer here 
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than for all other experimental maxima shown in Fig. 2. Except for asymmetry 
the agreement between experimental and calculated results is quite good con- 
sidering the many sources of disturbance in experiments on microwave 
diffraction. 

The experimental curves for the hard rubber rod in Figs. 3 and 4, if corrected 
for asymmetry and intensity by factors obtained from Fig. 2, should represent 
closely the true patterns of the rod in the chosen planes. The noticeable peak 
in the center of the shadow, which falls off in intensity with increasing x, is an 
effect of transmission and convergence by the cylindrical semitransparent rod. 
The rapid rise of intensity of the first peak outside the shadow in the lower six 
curves and the falling off in height of the next two in the same curves with 
increasing x are other distinguishing features of the patterns of the dielectric 
rod. Except for decreased peak intensities because of less scattering from the 
surface of the dielectric cylinder, the patterns in the planes x = 0 for it and 
for the conducting cylinder, Fig. 2, are closely similar. 

Diffraction patterns of the above cylinders and of a lucite rod (with negli- 
gible loss) are being studied with incident radiation polarized perpendicular 
to the cylinder axis. The diffraction field of the latter rod will also be studied 
using parallel polarization. Comparison with theory will be somewhat less 
difficult in the case of the lucite rod than in that of the hard rubber rod, which 
introduces considerable loss into the problem as an added complication. 
Results of these experiments will be reported at a later date. 


ACKNOWLEDGMENTS 


The authors wish to thank Mr. Bruce Chidley for making measurements of 
the dielectric constant and loss tangent of the hard rubber of the rod used in 
these experiments. One of us, A. B. McLay, thanks the National Research 
Council of Canada for award of a grant-in-aid for 1953-54 in support of 
continuation of this program of research. 


REFERENCES 


1. British Association Mathematical Tables, Vol. X., Bessel functions, Part II. Cambridge 
University Press, London. 1952. 

. Dakin, T. W. and Works, C. N. J. Appl. Phys. 18: 789. 1947. 

. Konis, R. D. J Appl. Phys. 23: 249. 1952. 

. Montcomery, C. G. Editor. Technique of microwave measurements. McGraw-Hill Book 
Company, Inc., New York. 1947. p. 57. 

Montcomery, C. G., Editor. Technique - phere measurements. McGraw-Hill Book 
Company, Inc., New York. 1947. p 

. Morse, P. M. Vibration and sound. MeGraw-Hill Book Company, Inc., New York. 1936. 


on m Gobo 


uo 


p. 
. STRATTON, J. A. Electromagnetic theory. McGraw-Hill Book Company, Inc., New York. 
1941. Chap. 6. 





THE ATTENUATION OF SECOND SOUND IN LIQUID HELIUM II 
ABOVE 1°K. ! 


By K. R. ATKINS AND K. H. Hart? 


ABSTRACT 


The second sound was in the form of a pulsed continuous wave with a pulse 
length of 1 to 2 msec. and a carrier frequency of 10 or 20 kc./s. The change in 
amplitude of the pulse was measured as the distance between the transmitter 
and the receiver was varied. To avoid boundary effects, no propagation tube was 
used and allowance had to be made for the spreading of the second sound beam. 
The attenuation was found to increase with increasing second sound amplitude. 
The attenuation extrapolated to zero amplitude had a finite value which increased 
rapidly as the temperature was lowered towards 1°K. Its order of magnitude was 
too large to be explained by viscosity effects, but was in good agreement with a 
thermal conductivity effect predicted by Khalatnikov. 


1. INTRODUCTION 


Thermal waves can be propagated in all substances, but are usually damped 
down in a distance of the order of one wave length. Second sound in liquid 
helium II is a thermal wave, but it involves a special process in which the 
‘normal and superfluid components oscillate out of phase with one another. 
This motion of the two components through one another encounters surpris- 
ingly little hindrance and the attenuation of second sound is, in fact, so small 
that it is not readily measured. However, there are two reasons for believing 
that an investigation of the attenuation would yield results of considerable 
theoretical interest. The first reason is that second sound provides a hydro- 
dynamical situation which is not complicated by the presence of solid bound- 
aries and therefore enables us to study effects which are independent of the 
boundaries. The second reason is that theoretical investigations such as those 
of Khalatnikov (9) indicate that the attenuation is at least partly due to 
basic processes involving the behavior of the elementary excitations which 
exist in the liquid. A study of the attenuation might therefore provide infor- 
mation about the nature of these excitations. 

Peshkov (15) deduced the attenuation from measurements of the band 
width of a resonant cavity at frequencies between 100 c./s. and 1600 c./s. 
Pellam and Hanson (14) have made similar measurements on a resonant 
cavity at frequencies in the neighborhood of 100 c./s., using a Rayleigh disk 
as a detector. In both cases the results are probably very dependent upon 
boundary effects. Pellam (13) also investigated the decay in the amplitude 
of a d-c. pulse of second sound as it travelled through the liquid. He found 
that the attenuation was very small below 2°K. but increased rapidly as the 
\-point was approached. Osborne (12), using a similar technique, observed 
much smaller attenuations in the vicinity of the A-point. The disagreement 
may be explained by the higher frequencies in Pellam’s 150 usec. pulse, as 
compared with Osborne’s 600 usec. pulse, or by the shock wave effects dis- 

1Manuscript received January 21, 1954. 
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covered by Osborne (12). They were both concerned with propagation inside 
a tube and so boundary effects may have been present in their experiments. 

With these considerations in mind, we have attempted to measure the 
attenuation under simpler and more sharply defined conditions. In order to 
retain the advantages of the pulse technique while working at a definite 
frequency, a pulsed continuous wave was used. The carrier frequencies of 10 
and 20 kc./s. were in the highest frequency range that could conveniently 
be used, in order to emphasize any part of the attenuation which might 
depend on the square of the frequency. The pulse length lay between 1 and 
2 msec. and included sufficient oscillations to ensure that the spread of the 
spectrum was only a small fraction of the mean frequency. To avoid boundary 
effects, no propagation tube was used and the second sound travelled freely 
through the liquid. Special attention was given to the possibility that the 
attenuation might depend on the amplitude of the pulse. 

A short preliminary account of this work has already been published 
(3). A detailed account, including those technical points which cannot be 
discussed here, is given in the Ph.D. thesis of K. H. Hart (8). 

2. APPARATUS AND METHOD 
2.1 The Second Sound Transmitting and Receiving Elements 

The transmitter and receiver were made by winding bare wire in zig-zag 
fashion on a plexiglas frame, with a spacing of about 0.05 cm. between adjacent 
portions of the wire. The transmitter wire was 80 » diameter chromel A and 
the receiver wire was 40u diameter phosphor bronze. As the second sound 
wave lengths used ranged from 1 to 2 mm., the frame was carefully designed 
to ensure that the wire was confined to a plane to within 0.2 mm. 

As shown in Fig. 1, the transmitter T was mounted above the receiver R 
and could be moved vertically through a measured distances The receiver 
could be aligned parallel to the transmitter by means of three rods such as J 
which extended up through O-ring seals in the cryostat cap. No reflection 
from any part of the arrangement was allowed to overlap the pulse travelling 
directly from the transmitter to the receiver. The active area of the trans- 
mitter or receiver was a square of side 2 cm. and the internal diameter of the 
helium Dewar was 5.4 cm. For the short wave lengths used the spread of the 
second sound beam was therefore too small for the walls of the Dewar to 
affect the propagation. 


2.2 The Electronic Circuits 

A block schematic diagram of the electronics is shown in Fig. 2. The pulse 
repetition frequency generator 1 produced triggering pulses at the rate of 
10 per second to start the oscilloscope sweep and to operate the gating circuits 
in the transmitter and receiver amplifiers. In unit 2 these triggering pulses 
were delayed by 1 msec. and then initiated a d-c. pulse which opened the 
gated amplifier 3 and passed through oscillations from the audio oscillator 4. 
The resulting pulse of continuous wave was amplified in unit 5 and fed through 
a transformer to the second sound transmitter 6. The peak power at the 
transmitter could be varied from 0 to 20 watts and was measured on the 
vacuum tube voltmeter 7. 
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Fic. 1. The cryostat and its contents. 
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Fic. 2. Block schematic diagram of the electronic circuits. 
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The oscillating Joule heat in the transmitter wire produced a temperature 
oscillation at twice the frequency of the input voltage and this travelled as a 
second sound wave towards the receiver 8. A steady current of 5 ma. was passed 
through the receiver wire and, as phosphor bronze has a high temperature 
coefficient of resistance at these low temperatures, the voltage across the 
receiver oscillated in unison with the temperature oscillations of the incident 
second sound wave. A step-up transformer passed this oscillating voltage 
to the input of the preamplifier 9 and thence to the main amplifier 10. The 
pass-band filter 11A was tuned to the carrier frequency and was incorporated 
to improve the signal to noise ratio. 

The object of the experiment was to measure the change in the amplitude 
of the received signal as the distance between the transmitter and the receiver 
was varied, but, as the attenuation was so small, this involved measuring 
changes in pulse amplitude with an accuracy of 0.1 db. or better. This was 
done by varying the attenuator 12 until the signal was equal in amplitude to 
a fixed comparison pulse derived from the input pulse to the transmitter, thus 
eliminating errors due to fluctuations of this input pulse. The attenuator 
was a carbon resistor network designed and constructed by Eckl (5). It was 
calibrated and found to be accurate to better than 0.05 db. in the frequency 
range from 1 kc./s. to 60 kc./s. 

The comparison pulse was taken from the transmitter input and was passed 
through the frequency doubler 13 to give it the same frequency as the signal, 
and then through the pass-band filter 11B which was as nearly as possible 
identical with 11A. The two pulses to be compared passed through the iden- 
tical channels 16A and 16B where they were amplified and clipped. A time- 
movable gating pulse from units 14 and 15 opened the signal channel 16A 
only during the passage of the desired signal, thus eliminating pick-up and 
unwanted echoes. The final comparison was made in a balanced integrating 
circuit similar to one used by Eckl (5). Each pulse was fed on to the plate of 
a 15 wF. condenser in parallel with a 1 megohm resistor, so that the condenser 
discharged by an amount representing the integrated effect of about 150 
pulses. The galvanometer G was connected between the two condensers and 
indicated the voltage difference between them. The procedure was therefore 
to adjust the attenuator 12 until the galvanometer was restored to its zero 
position. 


3. PRELIMINARY CONSIDERATIONS 


3.1 The Amplitude of the Received Signal 

Direct measurements of the signal amplitude were made by monitoring 
at a point immediately following the pass-band filter 11A, and measuring 
the height of the deflection on the screen of the cathode-ray oscilloscope 17. 
The gain of the amplifier being known, it was then possible to deduce the 
voltage swing across the phosphor bronze receiver. The temperature depend- 
ence of the resistance of the receiver wire had previously been measured and 
so the temperature oscillation, AT,, of the wire could be deduced. It is im- 
portant to consider to what extent this was different from the temperature 
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oscillation, AT;, in the liquid near the receiver. The equivalent circuit of the 
receiver has been given by Osborne (12) and is reproduced in Fig. 3. In our 
case the thermal impedance of the receiver wire, 1/C,w, was — 30 deg. cm.? 
watt—! at 20 kc./s., while the Kapitza boundary resistance, R,, at 1.2°K. is 


— dite Rt Transmission Line Rr 
To surface Y _1 of Characteristic 
of liquid Impedance Y 


Transmitter Liquid Helium Receiver 
Fic. 3. The equivalent circuit of the second sound apparatus. 


— 30 deg. cm.? watt~! (16), so the temperature discontinuity at the surface 
of the wire was probably comparable with the temperature oscillation in the 
wire. Another consideration is that the receiver wire deviated slightly from a 
plane and the second sound disturbances incident on different portions of it 
had slightly differing phases. Both these effects would be more important at 
higher frequencies and it is quite possible that, at 20 kc./s., the temperature 
oscillation in the liquid may have exceeded that in the wire by a factor of 
2 or more. 


3.2 Beam Spreading in the Liquid 

The temperature wave form in the transmitter wire was a consequence 
of the Joule heat developed by the input voltage and therefore consisted of 
a d-c. pulse superimposed on a pulse of continuous wave. The d-c. pulse 
component was composed principally of frequencies less than 1 kc./s., cor- 
responding to second sound wave lengths of several centimeters. This com- 
ponent therefore spread sideways very rapidly and was appreciably weakened 
at a distance of 2 or 3 cm. from the transmitter. Because of the nonlinearity 
of some of the attenuating mechanisms (Section 4.1), it was desirable not to 
have a steady heat flow present. As the attenuation measurements were all 
made at distances greater than 4 cm. from the transmitter, the d-c. component 
of the transmitter input could be ignored. The power dissipated by the steady 
current in the receiver and the natural heat leak into the cryostat were both 
small compared with the peak power of the second sound pulse. 

The carrier frequency of the transmitter input corresponded to a wave 
length in the range 1 to 2 mm., which was small compared with the width of 
the transmitter, and so this part of the second sound was sharply beamed. 
Most of the power was therefore incident on the receiver, but the small 


_sideways spread produced an apparent attenuation which was not negligible 


compared with the true attenuation. Apart from the fact that we were dealing 
with the theoretically difficult transition region between Fresnel and Fraun- 
hofer diffraction, departures of the transmitter wire from a plane rendered 
pointless any attempt to calculate the spread of the beam. The correction 
was therefore estimated experimentally in the manner described in the next 
section. 
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4. EXPERIMENTAL RESULTS 


4.1. Below 1.6°K. 

At a fixed frequency, temperature, and power input, the attenuator settings 
were recorded as the transmitter-receiver separation was decreased from 
9.5 cm. to 4.5 cm. The separation was always returned to its initial value 
to check that there had been no drift in the electronic circuits. This was 
repeated at several values of the peak power input P and curves of the type 
shown in Fig. 4 were obtained. The transmitter-receiver separation was then 






settings (10 log,, Alp) 
nN 
eee a eee 


I— 


Attenuator 


Transmitter - receiver separation - cm. 


Fic. 4. Variation of signal amplitude with transmitter-receiver separation. P = peak 
power input in watts. Frequency = 20 kc./s. Temperature = 1.20°K. 


set at 4.5 cm. and measurements were made of signal amplitude versus peak 
power input. For these measurements the signal from the pass-band filter 
was fed directly into the cathode-ray oscilloscope and the attenuator adjusted 
to produce a standard deflection of the trace. This procedure enabled the 
temperature amplitude in the receiver wire to be deduced for any power 
input and any transmitter—receiver separation in the range covered. 

A careful examination of the plots of attenuator setting against transmitter- 
receiver separation revealed that they were not quite straight lines. The 
curvature was too large to be explained by the variation of attenuation 
with amplitude and was probably a consequence of the superposition of beam 
spreading on top of the true attenuation. The total attenuation was taken to 
be a mean value over the whole range of separations and was derived by 
drawing a smooth curve through the points and taking the slope of the chord 
joining the two ends of the curve. The attenuation obtained in this way will 
be called the ‘‘total attenuation”’, a, and will be considered to correspond 
to the signal amplitude at the middle of the separation range. 
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It is obvious from Fig. 4 that the total attenuation increases with increasing 
power input and in Fig. 5, ais plotted against AT,, the zero to peak temperature 
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Fic. 5. Dependence of attenuation on amplitude at 20 ke./s. 


oscillation in the receiver wire. The vertical line drawn through each point 
shows the maximum uncertainty introduced by the slight curvature of the 
plots of the type shown in Fig. 4. Fig. 5 includes only results at 20 kc./s., 
but the effect was equally marked at 10 kc./s. To investigate whether the 
dependence of attenuation on amplitude was caused by some peculiarity in 
the electronic circuits, the power input was kept constant but the signal to 
noise ratio was varied by varying the steady current through the receiver 
wire. This left conditions in the liquid unaltered and involved no change in 
the nature of the electronic circuits. In this way the signal to noise ratio was 
varied from 1.5 to 10 without changing the apparent attenuation by more 
than + 3%, which was entirely random error. 
The total attenuation a can be divided into three parts: 


[1] a= a,+aot+ ay. 


aq is the “‘amplitude-dependent attenuation’’ caused by nonlinear dissipative 
effects in the liquid; ap is the ‘“‘zero-amplitude attenuation” which is analogous 
to an ordinary first sound attenuation in a normal liquid; a, is the ‘geometrical 
attenuation” caused by the spreading of the second sound beam. Extrapolation 
of curves similar to those in Fig. 5 to zero amplitude gives a, + ao, which is 
plotted against temperature in Fig. 6 for the 20 kc./s. results. The vertical 
line associated with each point represents an estimate of the uncertainty 
involved in extrapolating to zero amplitude. a, was estimated by assuming 
that ao is small compared with the experimental error between 1.5°K. and 
1.6°K. This is consistent with the theoretical predictions (4, 9, 10) and previous 
experimental investigations (12, 13). The value of a, obtained by smoothing 
the results between 1.5°K. and 1.6°K. was then taken to apply to the whole 
temperature range 1.1°K. to 1.6°K., since the velocity of second sound changes 
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Fic. 6. ag + ao at 20 kc./s. 


----- contribution from first viscosity ™.. ————— Khalatnikov’s prediction. _ 
(The two theoretical curves have beer moved upwards through a distance a, in order to fit 
at 1.6°K.) 


by only 7% in this range and the consequent small change in wave length 
should make no significant difference to the diffraction effects. The values 
of a, which had to be assumed are given in Table I. 
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----- contribution from first viscosity ,. ————— Khalatnikov’s prediction. 
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TABLE I 
THE GEOMETRICAL ATTENUATION 


Frequency, | Wave length at 1.5°K., 9, 
kc. /s. mm. cm,~! 
10 2 0.035 
20 1 0.018 


Fig. 7 is a plot of ao/w? where w is, as usual, 2x times the frequency. This 
graph shows that the zero-amplitude attenuation varies as the square of the 
frequency within the accuracy of our measurements. 


4.2. Near the \-Point 

Above 2°K. the signal to noise ratio was so small that it was possible to 
obtain results only at 10 kc./s. The few results obtained are given in Table 
II. At 2.05°K. there is an indication that the attenuation still depends on 


TABLE II 


ATTENUATION MEASUREMENTS ABOVE 2°K. 
FREQUENCY = 10 kc./s. 























| 
7, Wave length,} Power input, | AT,, a, Qo, ao + aa, 
*K mm. watts °K cm.~} cm.7! cm,.~! 
2.05 1.49 | 12.2 | 1.2 X 10-4 0.037 0.026 0.011 
2.05 1.49 20 1.7 10% 0.041 0.026 0.015 
| 0.95 X 107* 0.074 0.021 0.053 


2.12 1.16 | 20 





amplitude. The geometrical attenuation was estimated by assuming that it 
depended only on wave length and making a linear interpolation between 
the 10 kc./s. and 20 kc./s. values of Table I. There is a definite suggestion 
that the zero-amplitude attenuation is increasing as the temperature increases 
towards the A-point. Our values agree with those of Osborne (12) in being 
much smaller than Pellam’s (13) values, even though our frequency was 
greater than the frequencies in Pellam’s d-c. pulse. As the frequencies in 
the spectrum of Osborne’s pulse lay around 1 kc./s., it is satisfactory that 
his results lie below ours. 


5. THEORETICAL DISCUSSION 


5.1. The Amplitude-dependent Attenuation 

Although our results establish a variation of attenuation with amplitude, 
it is not easy to relate them to the various theories which have predicted this 
effect. The difficulty is that the experiments measure the temperature oscil- 
lation in the receiver wire, whereas the theories are concerned with the temper- 
ature oscillation in the liquid, and the factor of proportionality between the 
two varies in an unknown manner with temperature and frequency, as was 
already explained in Section 3.1. We are therefore restricted to order of 
magnitude arguments. 
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Three mechanisms which might give rise to an amplitude-dependent attenua- 
tion have already been postulated. They are: conversion to first sound, shock 
wave effects, and the Gorter—Mellink force of mutual friction. Breakdown of 
second sound into first sound has been observed by Pellam (13), but the 
breakdown was confined to a region within 0.25 mm. of the transmitter and 
it has not yet been established that the effect occurs in the free liquid. Pellam 
found that the breakdown set in above a critical power flow of 10-? watt 
cm.~*, which is certainly less than the power flows used in our experiments. 
Until the mechanism causing the conversion is understood no further discussion 
can be given; the relevant mechanism might, in fact, be one or both of those 
which will now be discussed. 

Shock waves of second sound were discovered by Osborne (12). They occur 
when the velocity of the normal component, v, becomes a sufficiently large 
fraction of the velocity of propagation, u2. In our experiments v, was never 
more than 20 cm. sec.~!, which is only about 1% of u2, but it does not follow 
that shock wave effects were unimportant. If it is permissible to apply the 
simple argument used by Osborne to explain the distortion of his pulses, 
then the peaks of our sinusoidal oscillation must be imagined to travel faster 
than the troughs and this would presumably produce a wave form of the 
saw-tooth type. A simple calculation shows that, with v, — 16 cm. sec.—!, the 
peaks would have advanced on the troughs by an appreciable fraction of a 
wave length in our experiments. A saw-tooth wave form was in fact observed at 
power inputs of about 20 watts, but we could not be quite certain that this 
was not due to distortion in the electronic circuits. Under these conditions, 
there was also an indication that the velocity of propagation had increased 
slightly, suggesting that shock fronts had indeed been formed, and at the 
same time a hiss at the pulse repetition frequency could be clearly heard 
coming from the cryostat. Apart from the attenuation known to exist in a sharp 
shock front, our experiment was particularly sensitive to pulse distortion, 
because the pass-band filter 11A passed only the fundamental frequency. 
These simple arguments suggest that the shock wave attenuation should 
depend on the frequency, inasmuch as it is easier for the peaks to overtake 
the troughs when the wave length is shorter. A thorough mathematical 
treatment is desirable, however. 

In order to explain the nature of thermal conduction in liquid helium through 
capillaries, Gorter and Mellink (6) postulated a mutual force of friction be- 
tween the normal and superfluid components having the form 


[2] Fy = Apspn(vs o Un)? 


with the constant A ~ 50. They pointed out that this should give rise to a 
nonexponential decay of second sound. If the attenuation per wave length 
is small, the work done against this mutual friction per cycle can be calculated 
and equated to the decrease in the total energy of the wave, giving an 
attenuation 


[3] a= 3 Apus (26 \'ars 
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where AT, is the zero to peak temperature swing in the liquid. The best case 
to discuss is the lowest frequency, 10 kc./s., and the highest temperature, 
1.6°K., because then the temperature oscillation in the receiver wire differed 
least from that in the liquid. Under these conditions, there was an attenuation 
of 0.018 cm.~! when AT, was 2.5 X 10-4 deg. K. Inserting this value for AT 
in equation 3 and taking A = 50, we deduce that ag = 0.006 cm.—'. The 
difference between the calculated and experimental a, has no significance in 
view of the possible difference between AT; and AT,, the uncertainty in the 
value of A, and the probability that F,, does not have the simple functional 
form of equation 2. Our results are therefore in order of magnitude agreement 
with the Gorter—Mellink theory. The accuracy is not good enough to establish 
that a, varies as AT®. 


5.2. The Zero-amplitude Attenuation 
Khalatnikov (9) has written the zero-amplitude attenuation in the form 





‘ ; 
[4] a= 2 [etd m+ fe — ob + 50) + o'r} +3] 
where 7, is the coefficient of first viscosity of the normal component, £j, 
fo, ¢3, and ¢4 are coefficients of second viscosity (actually ¢; = £4), x is analogous 
to the thermal conductivity of a normal liquid, and C is the total specific 
heat. Dingle (4), Kronig and Thellung (10), and Kronig, Thellung, and 
Woldringh (11) have obtained similar expressions, but Khalatnikov alone 
has attempted to give a physical significance to the coefficients £1, fo, £3, £4, 
and x and to calculate their magnitude. The dotted curves in Figs. 6 and 7 
show the contribution from 7,, using the experimental values of », obtained 
in oscillating disk experiments (1, 7). The observed attenuations are obviously 
significantly larger than this. The full curves show the complete expression 
in equation 4, as calculated by Khalatnikov, and the agreement with experi- 
ment is seen to be very satisfactory. The main contribution is from the thermal 
conductivity x, but the contributions from the first and second coefficients 
of viscosity are not negligible. The order of magnitude of x is 10-* watt 
cm.~! deg.—! near 1.2°K. 

According to Khalatnikov, the thermal conductivity x can be separated 
into two parts due to rotons and phonons separately: 
[5] X = Xr + Xpn- 
The two parts are approximately equal near the A-point, but xp, becomes 
progressively larger than x, as the temperature is lowered and is completely 
dominant below 1.4°K. A rearrangement of the formulae given by Khalatnikov 
. suggests that the viscosity and thermal conductivity may be expressed in a 
simplified form above 0.9°K. (2). For the rotons, 


Tv ~ 
[6] a. 10 Prvr Ar; 


[7] Xr = Br : purr, C.. 
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Here p, is the contribution of the rotons to the density of the normal com- 
ponent, 3, is the average velocity of the rotons, A, is the mean free path 
determined by roton-roton collisions, and C, is the roton specific heat. 8, is 
a factor approximately equal to 1 and varying only slightly with temperature. 
Similarly, for the phonons, 


[8] Npn = & Ppn C Ap; 


[9] Xph = Bon 3a pe ph Con 


where c is the velocity of first sound, @ is a constant not very different from 
unity, and 6, is also very nearly unity and varies very little with temperature. 
Although equations 6 to 9 are probably formally correct and help to demon- 
strate the nature of x above 0.9°K., it may be misleading to attach too much 
significance to them. Khalatnikov has emphasized that his x is analogous to 
a thermal conductivity, but is not quite the same thing as the thermal con- 
ductivity of a normal liquid. For example, Khalatnikov’s x for an ideal 
phonon gas is zero, whereas equation 9 suggests a very large value with Ap, 
limited only by the dimensions of the container. However, as far as the 
experiments on the attenuation of second sound are concerned, equations 
6 to 9 are adequate to explain the results so far obtained, and it must be left 
to further investigations to resolve the more subtle points. 


We should like to thank Professors J. M. Anderson and R. W. McKay for 
their assistance with the electronic circuits. During the later stages of the 
research, one of us (K. H. H.) was the holder of a scholarship from the 
Research Council of Ontario. 
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ELECTRON CAPTURE CROSS SECTIONS! 


By Harry SCHIFF? 


ABSTRACT 


Electron capture cross sections are calculated in the first Born approximation 
for alpha particles passing through hydrogen and singly charged helium ions 
passing through helium using the complete interaction Hamiltonian. Estimations 
of captures into excited states are made with the help of the partial cross sections 
obtained in a simple closed form ‘using only the (incident i ion) )-(electron) inter- 
action. The results indicate that the first Born approximation, using the com- 
plete interaction, is quite adequate in the velocity range given by e?/hu S 3. 
An impact parameter calculation for protons in hydrogen shows that the cross 
section obtained using only the (incident ion)-(electron) interaction gives 
unphysical results at low velocities (e?/fiv ~ 1), and that most of the contribution 
to the cross section arises from impact parameters p > do, where dp is the radius 
of the first Bohr orbit in hydrogen. 


INTRODUCTION 


In a previous paper (6) (referred to as I in the text), the electron capture 
cross section for protons in hydrogen was calculated in the first Born approxi- 
mation using the complete perturbation Hamiltonian. The results obtained 
were in good agreement with experiment for proton energies as low as 25 kev. 
(e?/hv = 1). Using the same methods, we now calculate the electron capture 
cross sections for alpha particles in hydrogen and, in a more approximate 
way, for singly ionized helium ions in helium. 

In Section 1 the results of Brinkman and Kramers (2) (referred to as BK 
in the text) are extended in a closed and more general form for purposes of 
comparison with the present calculations and as an aid in approximating 
cross sections for captures into excited states. In Section 2 the ground state 
capture cross section is obtained for the collision of a hydrogenic ion of charge 
ze with a hydrogenic atom of charge Ze, as well as the matrix elements for 
capture into the first excited state. In Section 3 the ground state capture 
cross section for protons in hydrogen is converted into an impact parameter 
form and is compared with the corresponding result of BK. Section 4 deals 
with the capture cross section for alpha particles in hydrogen. In Section 5 
an estimation is made of the capture cross section for singly charged helium 
ions in helium. Section 6 contains a discussion of results. 


SECTION 1 


Brinkman and Kramers calculated the ground state electron capture 
cross section for a hydrogenic ion of charge ze and a hydrogenic atom of 
charge Ze using only the (incident particle)—(electron) interaction as the 
‘perturbation causing the transfer. Following BK, the matrix element for the 
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capture into a final state f can be written in the center of mass system of the 
three particles as (see I, Eqn. (12)) 


[1.1] Hyon = — fexp (— 7B.r) ¢;(r) i exp (iA.g) ¢0(e) de dr, 


where ¢» is the ground state wave function for the electron in atom Z, @¢, is 
the wave function for the electron in the state f of atom 2, p is the distance 
of the electron from the nucleus, 7 is the distance of the electron from the 
incident particle, and 


cM bM 


AS attm *~ EB BOB itm © 
where, in terms of the initial and final relative velocities v and v’, 


hk = iV, hk’ = uyV’, 


uw, and uy being the reduced masses of the system before and after capture; 
cM and 6M are the masses of the nucleus and incident particle respectively, 
in terms of the proton mass M. 

For capture into a state characterized by the principal quantum number 
n, it is convenient to express the r dependent terms in the integrand of [1.1] 
in the parabolic coordinates, 

£ = r(1 — cos 8), 

n = r(1 + cos 6), 

% = ¢. 
With the choice of B as the polar axis, the dependence on the magnetic quantum 
number m is eliminated, and it is found that 


[1.2] a mcml = EE a ? 
: n1,ns(BK) (a + A*)* (a, + m2 + 1)(8°+ B’)*’ 


where 1, m2 are the parabolic quantum numbers and a = Z/a, B = 
2/do (my + m2 + 1), ao being the radius of the first Bohr orbit in hydrogen. For 
the particular choice of axes made, m; + m2 + 1 =, so that Eqn. [1.2] is 
independent of m; and m2, and depends only on m. Thus for the n™ state 
capture, 

25 10 54 2 2 

2° (aB)’e” ras 
1. Haw |? = Twn = 
[1.3] \Hncsx) | Z| .ns(BK)| (je +A?) (e+ BP 
where now, 8 = 2/ndo. With use of the conservation of energy equation 
(obtained from I) 
[1.4] a? + A? = pg? + B?, 
Eqn. [1.3] can be written as 


2510 542 4 
n2 (Zz/n)’e xa 
[1.3’] Hacer |” = 7 ap a 





Eqn. [1.3’] shows that the main contribution to the cross section comes from 
angles such that ao?A* < Z*. If ao?A? is expanded in a series in m/M one 
obtains 
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2 2 2 2\2 2 
42 2 eae Sey T, 
[1.5a] aA = 4 [(: ee +s rg ap ; 
also, 
[1.5b] a;B? =~ lise apiey » (ae )'] 
ie S\et+e/ J" 


where s = v/v, Yo being the characteristic velocity of an electron in the first 
Bohr orbit of hydrogen, and g = M@/m, where @ is the scattering angle in 
the center of mass system. After integration of [1.3’] over angle, the cross 
section for capture into the m™ state becomes 


[1.6] _ 2" (22/n)* nar's* 
= On(BK) = 5[s° “ (Z =e z/ny) [s+ (Z+2/n))° 


This result may also be written as 
[1.6’] on(BK) = n'o1(BK) (Z, z/n), 


that is, in terms of the ground state capture by a hydrogenic ion of charge 
z/n from a hydrogenic atom of charge Z. It is this relation which BK con- 
sidered as being approximately true only for high and low energies, whereas 
we have shown it to hold for all energies. For future reference we note that 


for s>Z+2/n, 


1 
[1.7] On(BK) ~ 53 F1(BK) 
for fixed Z, z. 
SECTION 2 


The complete matrix element (using both terms of the perturbation) for 

ground state capture is given by, 
i 2Ze ze 

21] Hi= fexp (— Br) a(t) ! ss at} exp (iA.g) $o(o) de ar. 
For this case the use of parabolic coordinates does not lead to any obvious 
simplification. The integration involving the (incident ion)—(electron) inter- 
action (—ze?/r) is straightforward. For the (incident ion)—(nucleus) integral 
it is convenient to express 1/|r ol, ¢o, and ¢; as Fourier transform integrals. 
Then one obtains the single k-space integral, 


io) at fo dk 
7+ (k — A)’]’ (6? + (k — B)’)* 





[2.2] Hy oue)1 = 


As shown in I, the integral in [2.2] can be evaluated with the aid of the 
auxiliary integral, 


J ‘  6x(1 — x) dx 
[2.3] ae ~ Jo jax-+ bd —x)) 


The complete matrix element can then be written: 
2 2ni(z, Z) —~ 1/2Z x (z, Z) INS 2, Za) 
ee a a 4. 4 


[2.4] Hy = 324(2Z)"e"ay : 
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where A =2Z?+a4,°A?, 
and : 
1 
x"(1 — x) dx 
& pa 
v(Z, = ae x ae? a ae ee 
A (2 Z) f — sa +} (s* “}- 3 a Z)x+ 3) / 


The integrals \5(z, Z) can be evaluated in a straightforward, though some- 
times tedious, way. 

If one puts v’ = v (the difference of the electronic binding energies being 
small compared to the energy of the incident ion), and expresses the dif- 
ferential solid angle in a small angle approximation, the cross section can be 
expressed as, 


[2.5] 1= (Ane'ag's')* f |H,|’ dA, 
th 


where ty = (4s?)—! [s? + (Z — 2/n)?] [s? + (Z + 2/n)?]. 
When the integration in [2.5] has been carried out, the ground state capture 
cross section becomes, 





[_V?, 2VW, W°+2VX , 2WX | 
— 7 2-2) VY sVW ae 

[2.6] o1 = 2°(2Z)' mao s + at} —— oe 242 + . 4h 
where V, W, X are the coefficients of A~*, A~?, A~! respectively in [2.4]. 
It can easily be shown that for z = Z = 1, [2.6] reduces to the result obtained 
in I for protons in hydrogen. 

For capture into the first excited state, the differential cross section is 
proportional to 

|H2|? = |Ho,|” + |Hop.|” + |Hop,|” + [Hop 

where the notation p,, p,, p, corresponds to the real hydrogenic p-state 
wave functions with azimuthal quantum numbers 1, —1, 0 respectively (8). 
The integrations for the matrix elements are carried out in the same way as 


for the ground state capture. The results are: 
2.29.4 
[2.7] Hp, ed (22Z)" xe*as E 52 Meo 2) oe 2/2Z 


2ni (2/2, Z) — ¢2°3(2/2, Z) — 1/2Z 
+ a See TE See ee 


4. M(e/2,Z) — iw 2"N5(@/2, Z) , 25(@/2, Z) — 9 2Mi(6/2, 2)| 


A A 





8(A3(2/2, Z) — 3(z/2, Z)) 
8(As(s/2, 2) — dals/2, 2) 
A 


12(A5(2/2, Z) — Ns(e/2, Z)) 4 15(A(E/2, Z) — M@/2, 2))| 
A’ ’ 


[2.8] Hop, = 12° (2Z)'2me'as A sin (A, B)| 





© A 


where, in a small angle approximation, one can write, 


; 2 oe, 
aed sin (A,B) = | s "42D ’ 
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[2.9] He,, = 0. 
(The results for Hz, and H2,, are obtained by a convenient choice of the 
coordinate axes for the k-space integration.) 
1 2(9)2 2 

[2.10] Hay, = i(22Z)""2meta’(A — 2°/4)74 = GNIS Beeld 2) 
4 3ri(2/2,Z) +3(s°+Z°—52"/4)A3(2/2,Z) —i(s°+Z°—2"/4)d3(2/2,Z) —1/2Z 

Soares ee nS FS ann eek earn hot 

A 
4 Aile/2, Z) + a(s* + Z* — 2°/2) N(G/2, Z) — 5(s* + Z* — 2°/4) ra(e/2, Z) 
ELST NE TE ire 


A 
‘ #r5(2/2, Z) + 33(s° + Z* — 2°/4)(AZ(2/2, Z) — dz(2/2, 2))| 
ceeemennieamminnaess whe snaccccks eee | 


It will be observed that the angular dependence of the matrix elements is 
contained only in the A’s. 

For the BK cross sections there exists the simple relation [1.6’] for the 
excited state captures in terms of the ground state capture. No such simple 
relationship holds when the complete interaction Hamiltonian is used. Since 
the complexity of the calculation increases rapidly with n, we shall only 
estimate the cross sections for captures into excited states with m > 2 (see 


Section 4). 
SECTION 3 


An enlightening comparison between the present calculations using the 
complete interaction and the calculations of BK can be made by expressing 
the total cross sections in each case in terms of an integral over all possible 
classical impact parameters. The integral is then just the probability that the 
capture takes place if the incident particle passes the atom at a given impact 
parameter. The possibility of conversion of our results into a meaningful 
impact parameter form depends on the fact that the nuclei are very heavy 
compared to the electron and their paths are essentially straight lines. Said 
in other words, the angular distribution is very strongly peaked in the forward 
direction; consequently partial waves of very high angular momentum must 
contribute to the cross section, and an approximate identification of a partial 
wave of given angular momentum with a classical path of a given impact 


parameter is possible. (See below.) 
The first step is to expand the matrix element [2.4] in spherical harmonics: 


[3.1] i we > (21 + 1) c:P;(cos 8), 
My i=0 


where k = yw,v/h is the reduced wave number. 
* The total cross section is then: 


oo 


(3.2] ¢= By (21 + 1) |e,|*. 


Since a partial wave of angular momentum /h corresponds approximately 
to particles with impact parameter p = //k, we can view the sum over / as 
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an integration of impact parameters ~, especially as a very large number of 
1 values contribute to the sum. Thus we can write the total cross section as: 


[3.3] o=2n J |\Crmup| b dp 


and ¢ metp|? is just the probability of capture at impact parameter p. 
To evaluate c; we use ss? [2.4] and [3.1] to obtain, 


[3.4] Cimtp = ah HP ;(cos 6) sin 6 dé. 


aa? 


The integrand is appreciable only at very small angles so that we can use 
the small angle approximation for P;(cos 8): 
P (cos 0) ~ Jo(l0) = Jo(py/ao), 


where y = (u,s/m)0, and we have used / = kp. From Eqn. [2.4] we see that 
we need integrals of the form: 


" Pi(cos 8) sin 6 dé <(2 57 © y Jo(py/ao) dy 


— woe AL. ae ynti 
0 A 


fl+s4t+ yy" 


with = 0, 1, 2. The last integral involves the explicit expression for A 
with z = Z = 1. This integral is readily evaluated (7) to give 


3 P (cos @) sin 6 dé 7 __x"Kn(x) 
s ys) 2"ni(1 + s°/4)”’ 


where K,, is the Hankel function of imaginary argument and of order n. 
x = (p/ao)(1 + s?/4)}. 
The complete expression for c;~;, is therefore: 
1 

13.5) crap = 29] fadca, 1) Kale) + MOD MEG) 4 a0. 1 sh a'Ka(s | 
The term in [3.5] independent of \% is the contribution from the (incident 
particle)—(electron) interaction and corresponds exactly to that obtained 
by BK in their direct impact parameter calculation, It can be easily shown 
that substitution of the above expression for c,-.) into equation [3.3] yields 
our previous result for the ground state capture cross section for protons in 
hydrogen (see Eqn. [2.6]). 

In Fig. 1 are plotted — c),, (indicated as P#) and the corresponding 
BK term as functions of p/a» for two values of the incident proton velocity, 
s? = land s? = 4. We note that whereas the BK curves are positive definite, 
the present curves change sign at about p/ay = 0.5 and p/ao = 0.3 at s? = 1 
and s? = 4 respectively, owing to the destructive interference of the two 
interaction terms. Also the BK curve for s? = 1 shows clearly that |c;|? > 1 
in the range 0 < p/ao < 2.25, having a value of about eight at p/ao = 0.25 
and decreasing to unity at p/a) = 2.25. This of course is physically meaningless 
since |c,|? represents the probability of capture and cannot be greater than 
unity. The BK curve for s? = 4 however is everywhere less than unity, having 
a maximum value of 1/2 at p/ao = 0. Considering the present curves in the 
same light, we see that they diverge (logarithmically as Ko(x)) for small 








| 
| 
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2 


Fig. 1. The spherical harmonic expansion coefficients, — ci-kp, of Hi in equation [3.1], 
indicated as P'/?, as a function of the impact parameter at two different velocities for protons 
in hydrogen. The broken curves refer to the calculations of Brinkman and Kramers and the 
solid curves are the present results. 


values of the impact parameter; the curve for s? = 1 is greater than unity 
for p/ao < 0.25 and has a maximum value of 0.55 at p/ay ~ 1.5, while the 
curve for s? = 4 is greater than unity for an appreciably smaller value of 
p/do than 0.25 and has a maximum value of 0.13 at p/ay) ~ 1. 

In order to examine the above observations in terms of the capture cross 
section, we have plotted, in Fig. 2, J = \e:|2p/ao, the integrand of [3.3] 


6 





— 10 I (S*=1) 

5 — 2501 (S*=4) 
Iex (S*=1) 
aa 25T,, (S*=4) 


2 
pe. 


Fic. 2. The probability of electron capture times the impact parameter (integrand of 
equation [3.3]) in units of ao as a function of the impact parameter at two different velocities 
for protons in hydrogen. The broken curves refer to the calculations of Brinkman and Kramers 
and the solid curves are the present results. 
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(divided by ao), for the same velocities as for the curves in Fig 1. The areas 
under the curves in Fig. 2 are proportional to the capture cross section. We 
observe here that for the BK curves most of the contributions to the cross 
section come from the impact parameter range 0 < p/ao < 3, the range 
decreasing with increasing velocity. As we have seen, it is essentially in this 
range that |c,|? is greater than unity for s? = 1, and thus the value of the 
BK cross section has no physical significance at this velocity. For the present 
curves we observe that in the region where P} diverges, the contribution 
to the cross section is very small and thus the cross section has physically 
plausible values for s? > 1. These results clearly show that the Born approxi- 
mation is more satisfactory towards lower energies when the complete inter- 
action is used. Finally we observe that the greatest contribution to the present 
cross section arises for values of p > do, the node in the curve moving towards 
smaller » values with increasing velocity. 


SECTION 4 
The results of Section 2 can be applied directly for the determination of 
the single electron capture cross section for alpha particles in hydrogen. The 
ground state capture cross section is obtained from [2.6] with z = 2, Z = 1. 
For the first excited state capture, the integration over angle of the modulus 
square of [2.7], [2.8], and [2.10] gives: 


2DE , E’+2DF , 2(EF + DG) 








: D® 
A aw tO aes as ee ae iu} 
[4.1] O25 < 7TQ0 sc 1z 7+ 6ts 6 5te° Ate 4 
F 4: 2EG | 2FG , G I 
2te" 





where ¢2 = (s? + 4)/4 and D, E, F, G are the coefficients of A-*, A-8, A-?, A-! 
respectively in [2.7]; 








[4.2] G2, = 2° mac | w +M+N)? ‘(1 = So og ae) 
ate L(L+2M) , 4(L ante i , item ¥) 
20t2* 122° 32,” 2te , 


where L, M, N are the coefficients of A~*, A~?, A~! respectively in [2.8]; 
st+4 ae HW , H(A +20) 
4 a Odlhl—— 
(H+ D* + 2HI , (H+ I)(H + I+ 2J) + 2HK 
5te” Ate" 
4 A+I+JI+K) —K’-2IK , (H+ JI4J34Ky -K 
= 2." 


,A+1+J+Ky ‘Tt | 


[4.3] oc, = 2!rae's *) (H +I+J+K)’ log 


o 





where /, J, J, K are the coefficients of A~*, A~*, A~*, A! respectively in [2.10]. 
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$* =(V/V,)* 


Fic. 3. Ratios of electron capture cross sections, as labeled, for alpha particles inhydrogen. 


In Fig. 3 we have plotted the curves of the ratios ¢2/0, and (o2/01)px in 
the energy range from 100 kev. to 1 Mev. (s? = 1 to s? = 10). We note that 
o2/o, varies from a value of 5 at s? = 1 toabout 1.1 at s? = 10. It can be shown 
that in the limit of high velocities o2/0; — 3, as is the case for (¢2/01)px. 
Thus, although the relative asymptotic contribution of a2 is small, at low 
velocities it is substantial, in agreement with the behavior to be expected 
from a resonance cross section (i.e., a cross section in which the electronic 
binding energy is the same before and after capture). 

It was shown in | that for protons in hydrogen there is a marked corre- 
spondence between the o2/o,; and (¢2/c1)sx curves, thus permitting an estima- 
tion of higher excited state capture cross sections in the terms of the BK cross 
section ratios. It is quite apparent in Fig. 3, however, that there is no longer 
the same correspondence between o2/o, and (¢2/c1)px at low velocities. 
Nevertheless, the higher excited states cross sections can be estimated in 
terms of the BK cross section ratios by the following considerations. 

For protons in hydrogen, since o; is the resonance cross section, the ratio 
o,/01 is always less than unity, whereas for alpha particles in hydrogen o2/o; 
varies from 5 to $; one might expect then a better correspondence between 
our ratios and those of BK if we compared o,/o2 and (¢,/e2)px for m > 2, 
since for alphas in hydrogen 2 is the resonance cross section and thus o,/o2 
for n > 2 is always less than unity. To show that this is so we have plotted 
(o2/03)px in Fig. 3. This curve varies from about 1.4 at s? = 1 to2.2 at s? = 10 
and we know that asymptotically it approaches (3/2)*. The position of the 
curve can be estimated as follows. First we expect that in the low velocity 
region, 1 < s? <4, where the resonance contribution of a2 is appreciable, 
o3 will be greater than o; since o3 is much closer to resonance than is 0}. 
Furthermore, ¢3 is certainly less than o2. Hence at low velocities, 1 < o2/a3 
<o2/o,and the lower bound is favored near s? = 1 since o; is appreciably greater 
than o;. An inspection now of Fig. 3 shows that with this restriction on 
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o2/o3, there is a rough correspondence with (¢2/¢3)sx. Also at high velocities 
we expect o2/0; — (3/2)%. Thus, to a fair approximation, o3 ~ o2(¢3/e2)sxK 
and it is reasonable to expect a similar behavior for the higher excited states. 
Approximately then, 
On & 02(0n/02) Bx = 02(03/02)BKx(on/o3)Bx- 

We have chosen the last expression for o, because, as an inspection of [1.6] 
shows, (¢,/03)nx ~~ (3/n)* for s? > 1 and » > 3, whereas (o,/o2)px varies 
appreciably with s? before attaining its asymptotic value (2/m)*. For the total 
cross section then we have 


Oe 0, +o. + 27 o2(03/o02) BK » 1/n° ’ 
or 
[4.4] oO, = a1 + ey + 2.1(es/ex)ax |} : 


In Fig. 4 are plotted o,, o2, and op) in the energy region from 100 kev. to 
1 Mev. Unfortunately, at present there are no experimental results for com- 


OC ClO cm?) 
i” 
- 
rai a 
oo 
¢ 
4 


go 2 3 4 $ 6 7 6 9 10 


Q-PARTICLE ENERGY (IOO KEV) 


Fic. 4. The total electron capture cross section, o-, (equation [4.4]) for alpha particles in 
hydrogen (solid curve). The cross sections for capture into the first excited state; present 
results — — —, results of Brinkman and Kramers - - - - - : 
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parison. It may be noted that o2gx) is greater than o2 by a factor of 8 at 100 
kev. and by a factor of 5'at 1 Mev. A comparison with the results of I shows 
that the cross section for capture of electrons by alphas in hydrogen is roughly 
six times that by protons of the same velocity in hydrogen, over the velocity 
range 1 < s? < 6. This factor of six can be explained in an approximate way 
as due to (1) the difference in charges of the alpha particle and the proton, 
(2) the increased number of final bound states in the helium atom per unit 
energy interval (see Reference (9)). 


SECTION 5 


We now give an estimate of the electron capture cross section for singly 
charged helium ions in helium. As an approximation for the wave function 
of each electron in the ground state of the helium atom, a hydrogenic wave 
function is used with effective nuclear charge Z = 27/16. As a basis for a 
further approximation, we recall the results of our impact parameter cal- 
culation in Section 3. It was shown there that most of the contribution to the 
cross section for protons in hydrogen arises from impact parameters p > ao; 
hence, insofar as the capture process is concerned we assume that the helium 
atom “‘sees”’ the incident helium ion with an effective charge of unity. 

Because of the conditions of resonance, the cross section for ground state 
capture will be by far the largest. Also for s? > 1 the collision time is short 
compared to the characteristic period of electronic motion. Consequently as 
a further approximation the effects of the uncaptured electron on the capture 
process are neglected. (This is analogous to the impulse approximation of 
nuclear theory (3).) Thus, to estimate the ground state capture cross section, 
it is only necessary to use the general one-electron result in [2.6] with 
z = Z = 27/16, and to correct for the magnitude of the interaction by dividing 
by (27/16)? since unit incident charge is assumed. To be somewhat more 
correct we should include an overlap integral factor for the two additional 
electrons because they both change states during the collision. This factor can 
be estimated to reduce the cross section by about six per cent, and will be 
neglected here. In terms of E = s?, the helium ion energy in units of 100 
kev., we find: 

3° rao $3.8443 . 0.72252 ie 
SEE + (27/8) Tl EF + EB + 291806 
9 3 ‘ & We 
- pS ae 





[5.1] oT = 


E + 1.9500 | 


ATs ae 7.6885 ; } 
a E | tan (2 2 a E + 2.6156 . 


Eqn. [5.1] is not multiplied by a factor of 2 for the two electrons in the struck 
helium atom because the Pauli principle allows the capture of only one of 
them into the ground state around the incident ion. 

For the excited states captures we assume that the relative contributions 
of the cross sections are roughly the same as for the asymptotic BK cross 
sections considered previously, that is, o, ~ (2/n*)o, (m > 2), the factor 2 


a 
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being due to the possibility of capturing either electron into an excited state. 
For the total cross section then, to a rough approximation, 


[5.2] o, = 01 (1 +2), 1/n’) = 1.401. 


In Fig. 5, we have plotted ¢, and two experimental points due to Snitzer (10). 
The agreement is remarkable and no doubt largely fortuitous, since for s? ~ 1 
the impulse approximation is inadequate. The agreement may not be as good 


100. 


O. (lOcn2) 


ae 2 3 4 5 6 7 
He ENERGY (IOOKEV) 


Fic. 5. Total estimated electron capture cross section (equation [5.1]) for singly charged 
helium ions in helium. The experimental points are those of Snitzer (10). 


for higher velocities, for then more contributions are to be expected from 
smaller impact parameters and the approximation to unit incident charge 
would deteriorate; our curve would then drop below the actual cross section 
values. That the theoretical curve is of the correct order of magnitude at 
s? ~ 1 can also be inferred from the experimental results (5) for s? < 1. 


SECTION 6 


We note that for the electron captures considered here as well as for protons 
in hydrogen, the charge of the incident particle was equal to or at most twice 
the nuclear atomic charge. Consequently, the determination of the resonance 
capture cross sections, which are dominant at low velocities, involved at most 
the calculation of captures into a first excited state (m = 2/Z). For cases, 
however, where the incident ion has a charge z > 2Z, the calculations would 
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be much more complicated since the determination of the resonance cross 
sections would then require the use of higher excited state wave functions. 
Furthermore, there would be more partial cross sections approximately at 
resonance contributing significantly to the total cross section. For such cases 
some approximation to the wave functions would probably be necessary. 

It is of interest to note that the validity of our expression [2.1], for the first 
Born approximation matrix element, has been justified by the work of Austern 
(1) and Gerjuoy (4) in their work on (d, p) stripping reactions. These results 
are sufficiently general so that the neutron capture in the (d, p) reaction is 
entirely equivalent, formally, to the electron capture problem. Examination 
of the integral equation obtained by these authors shows that the Born 
approximation using the complete interaction Hamiltonian is the first ap- 
proximation to the complete solution to the problem, and as was shown above 
(and in I), this seems to be quite adequate in the velocity range given by 
e*/hy 51. 

ACKNOWLEDGMENTS 


Acknowledgments are due to Professor J. D. Jackson for helpful discus- 
sions and assistance and to the National Research Council for financial 
aid in the form of a Fellowship. 


REFERENCES 


. AusTERN, N. Phys. Rev. 89: 318. 1953. 

. BRINKMAN, H.C. and Kramers, H. A. Proc. Acad. Sci. Amsterdam, 33: 973. 1930. 

CueEw, G. F. and Wick,G.C. Phys. Rev. 85: 636. 1952. 

. Gerjuoy, E. Phys. Rev. 91: 645. 1953. 

. HastepD, J. B. Proc. Roy. Soc. (London), A, 205: 421. 1951. 

; JACKSON, 4b. and ScuiFF, H. Phys. Rev. 89: 359. 1953. 

. MaGnus, and OBERHETTINGER, F. Formulas and theorems for the special functions 
of mathematical hysics. Chelsea Pub. Co., New York. 1949. 

. PauLinG, L. and Witeon, E Introduction to quantum mechanics. McGraw-Hill 


Book Company, Inc., New Voi and London. 1935. 
. ScuiFF,H. Ph.D. Thesis, McGill University, Montreal, Que. 1953. 
. SNITZER, E. Phys. Rev. 89: 1237. 1953. 


_ 





DEAD-TIME LOSSES IN MULTICHANNEL COINCIDENCE 
SYSTEMS! 


By J. S. KirKALpy 


ABSTRACT 


Correction formulae for the counting losses due to dead-times in multichannel 
coincidence systems are derived for conditions of significant interchannel 
correlation and for both continuous and pulsed sources of counts. The dead-times 
are assumed to be either ‘‘extending’’ (following all counts) or ‘“‘nonextending”’ 
(following only recorded counts). The corrections for the effects of a finite 
coincidence resolving-time are also evaluated. The approximations used are 
shown to be valid when the fractional losses are less than 40%, by comparison with 
the results of a previous experimental paper. The corrections are most likely 
to be required in the case where the source is intermittent, as from a pulsed 
accelerator, so that the instantaneous counting rates are very high during the 
pulse, and where the coincidence rate is of the same magnitude as the non- 
coincident rates, as in a well-shielded coincidence telescope. 


1. INTRODUCTION 


Westcott (5) and Westcott, Greenberg, and Kirkaldy (6) (hereafter referred 
to as I and II respectively) have studied the counting losses in single-channel 
and two-channel coincidence systems with a ‘‘nonextending” dead-time. 
In II an empirical formula was obtained for the losses in coincidence counting 
under conditions where the total counting rates in the respective channels 
do not greatly exceed the coincidence counting rate. In the present paper an 
approximate theoretical relation is obtained for this case, the counting rates 
all being constant (‘‘unpulsed’’), giving the fraction of coincidences recorded 
as a series expansion in ascending powers of the dead-time.* It will be shown 
that the retention of terms up to the second order in this series gives an 
accuracy of better than 1% in those cases where the fractional losses do not 
exceed 40%, this condition being usually met in practice. To this order the 
results are easily extended to systems of more than two channels. Corre- 
sponding solutions are also given for the alternative ‘‘extending’’ dead-time 
mechanism and for systems excited by a pulsed accelerator. 

We consider the dead-time 7) to be constant and to follow either all counts 
(‘‘extending”’ dead-time) or only those counts which record (“‘nonextending”’ 
dead-time) (4). A rectangular coincidence resolution curve of width 27, is 
assumed, 7, being the coincidence resolving-time. In cases of pulsed counting, 
the source pulse profile is assumed to be rectangular and of constant length -77 





2. TWO-CHANNEL COINCIDENCE COUNTING LOSSES (7, = 0) 


2.1. A ‘“‘Nonextending’’ Dead-time 
For a single channel with this mechanism operative and with an average 
rate of occurrence of events 7o, we can write the rate of recording counts, 7, 


1Manuscript received March 11, 1954. 
Contribution from the Radiation Laboratory, McGill University, Montreal, Canada. 


*More precisely, the series is an expansion in ascending powers of a physically dimensionless 
quantity, t.e., a product of a counting rate and the dead-time. 
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as ro multiplied by the fraction of the time that the channel is not dead, i.e., 
as ro(1 — rro). Solving, we obtain for the fraction of counts recorded (2) 


[1] at _, 
To ] ok ToTo 
If we consider two coincidence channels with rates a + ¢ and 6 + ¢ respec- 
tively, the ‘‘c’’ counts being coincident, we may approach the problem by 
evaluating the fraction of the time during which neither channel is dead, 
thence obtaining the fraction of coincidences recorded by analogy with [1]. 
Fig. 1 illustrates the types of occurrences that affect the result to the second 


c a a a ae ¢ 
(ne + «ee ee 

| l ' 

i | | 

' b b b 1 b | 


CHANNEL 2 J Fo a ee) Sa 
Fr 
CASE ! 2 3 4 5 6 " 


Fic. 1. Illustration of the contingencies involved in two-channel coincidence counting. 
Hatched sections represent overlap times. 


order in the dead-time (cf. footnote p. 406). We may therefore write the 
fraction of the time that neither channel is dead to the second order as 


[2] No/e = 1 — {Noto + pidto + pobto — pipsabro? + pip2(a + b)cro°/2} 


where p; = 1/{1 + (@ + c)ro} and ps = 1/{1 + (6+ c)ro} are the single- 
channel recording probabilities given by [1]. N2 is the rate of recorded coin- 
cidences (to be evaluated) and contributes an amount N27» to the fraction 
of the time either or both of the channels are dead (case 1 of Fig. 1). The 
second and third terms within the bracket include the dead-times due to all 
recorded ‘‘a” and ‘‘b’’ counts (cases 2, 3, 4, 5, 6, and 7) from which is sub- 
tracted the fourth term representing the fraction of the time during which 
“a” and ‘“d"’ counts overlap (cases 4 and 5). The fifth term represents the 
additional fraction of the time due to ‘‘c’’ counts recorded in one channel 
only (cases 6 and 7). Solving for N2/c and simplifying, we obtain for the 
fraction of coincidences recorded 

[2a] N2/c => pip2{ 1 ++ CTE = (a -- b)cro?/2}. 


To extend the above calculation to the third order we make use of elements 
of an alternative formulation of the problem (cf. § 5.1 of II). Following 
Fig. 2 we specify that a count (‘‘a” or ‘“‘c’’) has recorded at time ¢ = 0 in 
channel 1 and ask what is the conditional probability, P.(¢), that a count 
will record at time t (—t>) < ¢ <0) in channel 2. This is given by Feather’s 
relation (3) 


[3] P(t) =1- P.(t) Fo(t) dt 


t—To 
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Fic. 2. Illustration of the relations between the probabilities involved in coincidence 
counting. 


where F;(t) is the conditional instantaneous rate of occurrence of events 
in channel 2. The main problem is to determine F,(t) within the range 
—t) <¢ <0; this must be significantly less than its normal value b+ c¢ 
on account of the prior condition that a count has recorded at ¢ = 0 in channel 
1, since a “‘c’’ count can only occur in channel 2 in the interval to before 
t = 0 if its counterpart in channel 1 has been suppressed by some count 
recorded prior to ¢ = —7o9. On the average, to suppress a ‘‘c’’ count in the 
interval — ro <¢ <0, an “a” or “c’’ count must record in the interval 
— $31) <t< —7; the probability of this is approximately (a + c)ro/2 so 
that the conditional rate of occurrence of suppressed ‘‘c’’ events in channel 
1 in the interval —ry) < ¢ <0 is approximately c(a + c)ro/2. In evaluating 
F,(t) we drop the c?ro/2 term since, to the order of this calculation, a ‘‘c’’ count 
suppressed in channel 1 by another ‘‘c”’ is also suppressed in channel 2 and so 
can contribute nothing to the value of P2(t) in the interval. We therefore take 
F,(t) = 6 + caro/2 for —1t) < t < 0, F(t) = b+ cand P,(t) = p2 fort < —7. 
Differentiating [3] and solving (cf. § 5.1 of II), we obtain to the second order 
int 

[4] Pe(t) = po{1 + c(ro + t) — (@ + d)c(ro + #)?/2}, —m<t<0.* 

If the initial count taken at ¢ = 0 is a “‘c’’, the conditional probability of the 
counterpart in channel 2 recording is P2(—e) where e is an infinitesimal 
quantity. Since the probability of the initial event recording is pi, the expected 
fraction of coincidences recorded is 


N2/c = pipe {1 + cro — (a + b)cro?/2} 


which was the result obtained in [2a]. This formulation is very convenient 
for evaluating terms up to the second order and will therefore be used in 


*4 rigorous (and rather involved) calculation shows this relation to be correct throughout the 
range to the first order only, but as t/ro— O in the same range the relation becomes correct to the 
second order. 
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our subsequent calculations for systems with more than two channels. The 
third-order terms for a two-channel system (the only one considered to this 
accuracy) are most conveniently evaluated by extending our initial method 
(cf. equation [2]). Only the last two terms of this equation have to be modified. 
The fraction of the time that ‘‘a” and ‘d’”’ counts overlap (cases 4 and 5 
in Fig. 1) is to the third order 


0 
an f bP2(t) (t a. To) dt = abpipeto (1 — § cro) 


where P,(t) is taken to the first order in ¢ from [4]. A similar procedure applies 
to cases 6 and 7 and to four other contingencies which involve the overlap 
of three recorded counts and contribute to the third-order terms. The fraction 
of coincidences recorded is thus obtained to the third order in 7) as 


[5] N2/¢ = pipe {1 + ¢to — d(a + b)cre? + 4(a? + 4ab + b*)cro*}. 


Fig. 3 shows the second-order (equation [2a]) and third-order (equation 
[5]) forms plotted against tr) and compares them with the empirical equation 
and experimental points of II fora particular case. From this and similar 
comparison curves we conclude that the second-order form gives results 
accurate within 1% provided the losses are not greater than 40%, while the 
third order is accurate within 1% provided the losses are not greater than 60%. 
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Fic. 3. Fraction of coincidences recorded as a function of the dead-time. 
a = 184.2 cts./sec. b = 174.8 cts./sec. c = 265.6 cts. /sec. 
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2.2. An ‘‘Extending’’ Dead-time 

In a single channel with this counting mechanism, the probability of any 
count recording is the probability that no other count has occurred in the 
interval 7) prior to the count in question. Hence if the rate of occurrence of 
events is ro then the expected fraction recorded is 
[6] r/’o = e-Te7, 


For a two-channel coincidence system the analysis is similar. The probability 
that any “‘c’”’ count records in channel 1 is exp {—(a@ + c)ro} and the prob- 
ability that its counterpart in channel 2 also records is the probability that no 
6” count has occurred in the interval to prior to the count in question. A 
‘‘c”’ count is excluded from channel 2 in this interval by our prior condition 


in channel 1. The fraction recorded is thus 
[7] N2/c = e(atote) ro 


which is the same, to the first order, as our expression [2a] for a ‘‘nonextending”’ 
dead-time. 
3. MULTICHANNEL COINCIDENCE COUNTING LOSSES (7, = 0) 
3.1. A ‘‘Nonextending’’ Dead-time 
It is not difficult to generalize the two-channel results of § 2.1 to multi- 


channel systems. We consider in particular the “telescope” arrangements 
of Fig. 4. Fig. 4a shows an arrangement consisting of detectors which have 


D, D2 Ds 
Qa; G2 Qs Qh-1 On 


Fic. 4. Coincidence telescope arrangements; (a) an ‘‘ideal’’ m-channel telescope, (b) a 
“practical” three-channel telescope. D = detector; S = source; other symbols are counting 
rates. 


individual noncoincident rates a; (¢ = 1, 2,.. . ) (perhaps due to background) 
and an n-fold coincidence counting rate c. We call this an ‘‘ideal’’ system, 
since we assume there are no coincident events which excite less than n 
detectors. Hence the first » — 1 detectors must be perfectly thin or the 
energy of all the source particles must be sufficient to traverse all detectors. 
Fig. 40 is a “practical” three-channel system in which we allow for the events 
neglected above. Rates a, b, and d represent counts in one detector only. 
Rates c,; and care source (or possibly background) coincidence rates exciting 
the first two and all three detectors respectively, while rate cz represents 
coincidences between the last two detectors when the first is not excited. 

For the “ideal” system we evaluate the recording probability to the second 
order following the methods of § 2.1. We note that the probability that ‘‘c’’ 
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counts are recorded in both of the first ‘two channels is pipe {1 + cro 
— c(a; + d2)r¢?/2} (cf. equation 2a) and we require the conditional probability, 
P;(t), that these counts also record in channel 3. Noting that to the order 
of our calculation the conditional rate of occurrence of counts in channel 3 
is F3(t) = a3, —t) <t<0, and F;(t) =a3;-+c¢, t < —70, we can solve 
the relation 


Pit) =1- f Pe) Oat 


to obtain 
P;(—e) = ps{l + CTo — ca370?/2}. 


For each following channel in turn we can write 
P;(-€«) = pi {1 + ct) — ca;7 9/2} 


and the product of all these (¢ = 3,4, . . . m) and the initial two-channel 
recording probability yields the expected fraction of ‘‘c’’ events recorded in 


n channels, 
f 2 


n n 
[8] N,/c = it + (n —1)cre+ So Oe 2s “o a, aT] ais 
2 2 2 t= 
It seems reasonable to apply the same bound for less than 1% error as in the 
two-channel case, viz. N,/c > 0.6, since in this way we demand that the 
allowable counting rates be decreased as the number of channels is increased. 
For the ‘‘practical’’ system the probability of a coincidence being recorded 
in the first two channels (irrespective of whether the count is a ‘‘c;” or a ‘‘c’’) 
is given by [2a] as 


Pip2 {1 - (e + c¢)ro — (1 + c)(a +b+ C2)T0°/2} 


where the p’s are the appropriate single-channel recording probabilities. 
The conditional probability, when a coincidence is known to have recorded 
in the first two channels at ¢ = 0, that a count in channel 3 records at time 
t(—t) <t <0) is given by the integral equation above, where F;(t) is 
d+c+cefort < —7,and remains to be determined in the range —t) < t < 0. 
The rate d is unaffected by our prior condition in the latter interval, but a 
‘co’? may only occur (neglecting second-order contingencies) if suppressed 


in channel 2 by a “‘b” or ‘‘c;’’ which records in that channel prior to t = —79 
and a ‘‘c’? may only occur if suppressed by a ‘‘c;’’ which records in both 
channels 1 and 2 prior to ¢ = —7o. We therefore take 

F;(t) = d + 3e:cro + 3(b + €1)CaT0, —m<t<0O, 
and P;(t) = p3, t < —7 9, and solve to obtain P;(¢). Enumerating for the 
““c’’ coincidences only (i.e., taking ¢ = —e) we finally obtain for the expected 


fraction of triple coincidences recorded to the second order 
[9] N3’/c = pipops {1 + (cr + co + 2c)ro + (cr + C2 + 2c)er?/2 
—[(a + d)er + (6 + dco + (@ +b + d)c]re?/2} 


with the proviso for 1% accuracy that N3'/c > 0.6 as before. 
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3.2. An ‘‘Extending’’ Dead-time (rt. = 0) 

The evaluation of the fraction of coincidences recorded in multichannel 
systems is a straightforward generalization of the method outlined in § 2.2 
for two channels. For the ‘‘ideal”’ system 


[10] N,/¢ = exp — a + ct To 


and for the three-channel ‘‘practical”’ system 
[11] Ny! /c = e~atbtaterteste) re, 


To the first order these relations are equal to the corresponding ‘‘non- 
extending’’ dead-time results ({8] and [9]) but differ significantly in the higher 
order terms. 


4, PULSED COINCIDENCE COUNTING 


The calculations for the losses in one- and two-channel systems with a 
‘“‘nonextending”’ dead-time and excited by a pulsed source with a pulse- 
length Z have been carried out in I and II. We extend these methods below 
to multichannel systems. 

In general we may represent the number of counts recorded ina time J as 


TF 
J P(t) F(t) dt 


where P(t) is the instantaneous probability of a count recording and F(t) 
is the instantaneous rate of occurrence of events. In any of the coincidence 
systems previously considered, with a pulse-length FZ < ro, P(t) = exp(—xt?) 
is the probability that no count has occurred within the pulse prior to ¢ 
where x is the sum of all counting rates, coincident and noncoincident, ap- 
pearing in the system (e.g. for two channels x = a + 6 +c). Fora rectangular 
pulse F(t) is a constant (in the two-channel case F(t) = c), so on evaluating 
the integral and dividing by the expected number occurring during the pulse, 
J F(t), we obtain for the fraction recorded in m channels 
—17T 


l-e 
[12] on(T) = TGF) F<. 
For a two-channel system this gives 
— e atts 


(atb+cZ7 


and similar expressions may be obtained from [12] for the systems with more 
than two channels considered previously. These results apply independently 
of the dead-time mechanism considered, since at most one count can record 
per pulse. For 7 > 79 we can approximate, provided the losses are not too 
great, by taking P(t) for t < 7) as equal to exp(—xt) and for t > 79 as equal 
to the recording probability for continuous counting. This has been shown 
to be valid within reasonable bounds for a single channel (cf. Fig. 3 of I) and 
is valid within similar bounds for multichannel systems. Hence for any of the 


[13] ¢:(7) = (cf. II, §8.1) 
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coincidence systems previously considered, whether an ‘“‘extending’’ or 
“‘nonextending”’ dead-time is operative, we may perform the above integration 
to obtain the fraction of counts recorded 


[14] dn(J) = us +24 da(or) = Nat , i 


where N,/c is the appropriate continuous counting recording probability and 
¢n(To) is given by [12]. Fig. 5 shows a particular plot of ¢2 (given by [13] and 
[14] with N2/c from [2a]) against the pulse-length in units of the dead-time 
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Fic. 5. Fraction of coincidences recorded as a function of the pulse-length in units of the 
dead-time. a = 181.4 cts./sec. b = 188.0 cts./sec. c = 200.5 cts./sec. ro = 0.95 msec. 


(7/7) along with the corresponding experimental results from II. The agree- 
ment is seen to be quite good in the region of measurement. It appears from 
these and our other results that reasonably accurate corrections for dead-time 
losses may be obtained by the above formulae if we specify that N,/c > 0.6, 
a bound which is already required to assure the accuracy of our second-order 
approximations for N,/c in §§ 2.1 and 3.1. 


5. EFFECTS OF THE COINCIDENCE RESOLVING-TIME 


In II, it was shown that the effect of the coincidence resolving-time 7, 
on the rate of recording coincidences in a two-channel ‘‘nonextending” dead- 
time system can be correctly represented by a multiplying factor, viz. 


[15] N2(ro, te) = No(ro, 0){1 + (a + 6 + 2ab/c)r-} 


where N2(7o, 0) is given by equation [2a] or [5]. This relation gives good 
agreement with experiment provided 7, < 0.279 (cf. Fig. 13 in II for the 
experimental results). The ‘‘a” and ‘‘b”’ terms are due to one partner of a 


ere rears 


a RS Fr 
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coincident pair being suppressed by a noncoincident count which records 
in the interval 7, prior to the coincident event (a ‘“‘restored”’ coincidence) 
and the last term is due to an “‘a”’ and ‘‘b” count recording within an interval 
te (a ‘‘random” coincidence). A similar relation can be written for the ‘‘ideal’”’ 
multichannel system. For n > 2, first-order terms may only arise from 
coincidences recorded in any m — 1 channels and “restored” in the n*. 
Second and higher order terms are generally negligible for 7. < 0.279 and 
less than 40% losses, but an exception must be made for the term describing 
“‘random”’ coincidences between the ‘‘a;”’ counts. Although this term is of 
(n — 1)™ order in our expansion it may be numerically predominant when c 
is quite small compared with the noncoincident rates. To include this case 


we therefore write for the rate of oe coincidences 





n—l n 
[16] Ni(ro, 2) = Nn (ro, 0) + y agre+ a “71 a} 
i= 


where NV,,(r0, 0) is given by equation [8]. A simple derivation of the last term 
is given by Arley and Rander Buch (1). It is important to note that the cor- 
rection term in the above form is not a series expansion in the usual sense, 
since there is no restriction on the magnitude of the last term (cf. II, Fig. 
13c, which shows a particular case for » = 2 where the latter term has a 
value — 5 in agreement with the experiments). In the common case, where c 
is small relative to the other rates, this term alone is important and the 
requirement rt, < 0.27) may usually be relaxed to 7, < 0.579. Further relaxa- 
tion of the bound (a requirement which is unlikely in practice) introduces new 
possibilities, e.g., for 7 = 2 and r, > 0.579, two counts recorded in one channel 
may both lie within 7, of a count recorded in the other and the question arises 
whether in this event one or two coincidences record. If these higher order 
contingencies are negligible and c is small we may discard the bound altogether. 

For the “practical” three-channel system we can obtain by similar methods 


[17] Ns’ (ro, 7.) = Ns’ (r0,0)}1 A. (tetera ees ae 





with similar bounding conditions to the case above. The correction terms 
in equations [16] and [17] apply also for ‘‘extending” dead-time systems as 
well as for pulsed counting provided in the latter case, Z > ro, Nn/c > 0.6, 
and tr, < 0.179 (cf. II, § 9). 


“é 


6. DISCUSSION 

It must be pointed out that in many practical cases the coincident counting 
rate is much smaller than at least one of the noncoincident rates; the channels 
are then effectively independent so that our counting loss formulae reduce to 
simple products of the single-channel recording probabilities. Our complete 
formulae will probably only be required in the case where the source is inter- 
mittent, as from a pulsed accelerator, so that the instantaneous rates are 
very high during the pulse, and where the coincident rate is of the same magni- 
tude as the noncoincident rates, as in a well-shielded coincidence telescope 
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arrangement. In this case the formulae presented will give accurate corrections 
provided we make allowances if the dead-time mechanism is not the ideal one 
postulated and, in pulsed counting, if the pulse envelope is not rectangular. 
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SEARCH FOR A POSSIBLE ERROR IN THE 
MEASURED HALF LIFE OF Au’?! 


By R. E. BELL? AND L. YAFFE 


ABSTRACT 


The half life of Au’ has been measured using a quartz fiber electroscope, 
with extreme precautions against errors due to the presence of Au’ in the source. 
The result obtained, 2.699 + 0.003 days, shows that other recent determinations 
of the Au’ half life do not suffer an appreciable error from this cause. 


INTRODUCTION 


An accurate knowledge of the half life of Au'* is necessary when the cross 
section of the thermal neutron reaction Au!”(n, y)Au'® is to be used as a 
standard of comparison, as well as when Au’ is used as a comparison source 
in beta-ray spectroscopy. Several determinations of this half life have been 
made with good precision, as summarized in Table I. The statistically weighted 


TABLE I 


SUMMARY OF RECENT MEASUREMENTS OF THE HALF LIFE OF Au!% 
References to those in the upper sections of the table may be found in Hollander et al., Ref. (4) 























Duration of Half life 

Radiation detector measurement obtained 
Authors employed in half lives (days) 
Diemer and Groendijk Geiger counter 2.5 2.73 + 0.02 
Silver Ionization chamber ~12 2.69 + 0.01 
Saxon and Heller Geiger counter 5 2.69 + 0.02 
Sinclair and Holloway Geiger counter 4 2.72 +0.01 
Sinclair and Holloway Ionization chamber 3 2.74 + 0.015 
Cavanagh et al. Geiger counter 7.5 2.66 + 0.01 
Weighted mean of above values 2.698 
Probable error from quoted individual errors 0.003 
Probable error from consistency of results 0.008 
Lockett and Thomas (5) Ionization chamber 4 2.697 + 0.003 
Present work Ionization chamber 10 2.699 + 0.003 


mean of the values in Table I is 2.698 days, with probable error 0.003 days 
computed from the quoted errors of the individual entries, or 0.008 days 
computed from the scatter of the individual entries about the mean. Half 
life measurements on Au! are, however, subject to a possible error which is 
peculiar to this particular activity, as follows. Hill and his co-workers have 
shown that Au’ has an unusually high neutron capture cross section, accord- 
ing to their latest results about 3.8 X 104 barns (3, 7). Thus, for example, 
if gold is irradiated to saturation in a thermal neutron flux of 10'* cm.~ sec.“, 
there will be produced 12.8% as much Au’ as Au"’, by successive capture of 
neutrons. The reactions are 


1Manuscript received March 9, 1954. 

Contribution from Radiation Laboratory and Department of Chemistry, McGill University, 
Montreal, Quebec. Issued as A.E.C.L. No. 102. 
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Au" (n, y)Au™(n, y)Au™, 


The half life of Au’, 3.15 days (1), differs by only 17% from that of Au'®, 
so that a decay curve of the mixed source of Au’ and Au" is unlikely to be 
resolved into two components; instead, the initial slope of the decay curve 
will show a half life, in our example, about 2.2% too high. Owing to the 
slightly slower decay of Au’, the error in slope gradually increases with time 
after the end of the irradiation. The magnitude of the error in any particular 
case tends to be proportional to the neutron flux in which the Au'® source 
was prepared. 

(An additional way in which Au’ could be produced in an Au" source 
involves platinum impurity in the original gold sample. The reactions are 

B- 
> Au™. 
31 min. 


Pt*(n, 7) Pt} 





Consideration of the likely levels of platinum impurity in gold shows that this 
source of érror cannot be important.) 

We have measured the decay rate of Au" over an extended period, taking 
precautions to ensure that Au was not present in the source in quantities 
large enough to cause an appreciable error in the Au half life. 


EXPERIMENTAL 


Two similar pieces of spectroscopically pure gold foil (platinum content 
less than 1 p.p.m.) were irradiated in the thermal column of the Brookhaven 
reactor for 24 hr. at a mean thermal neutron flux of 7 X 10° cm. sec.“. 
One of the foils was wrapped in 1 mm. cadmium, and had an activity of 
1/3000 of that of the main sample. Thus the activity of the main sample 
was almost entirely due to thermal neutron capture. A small section of the 
main sample was cut off and checked for the presence of Au by searching 
for the prominent 158 kev. gamma ray (3, 7) of Au with a sodium iodide 
scintillation spectrometer. The main gamma radiation present was the 412 
kev. ray from Au’, A small peak was present in the pulse height spectrum 
at about 158 kev., which, if attributed entirely to Au, would indicate an 
admixture of about 1.3%. This peak, however, occurs at an energy corre- 
sponding to 180° Compton scattering of the 412 kev. gamma ray from the 
photomultiplier and crystal mounting. The peak was shown to be mainly 
due to backscattering by repetition of the scintillation spectrum with a 
4.2 gm./cm.? lead absorber shielding the counter. This absorber would reduce 
the intensity of a real 158 kev. gamma ray by a factor 1/550 relative to the 
intensity of the 412 kev. gamma ray, but in fact no relative reduction was 
detected in the peak in question, to an accuracy of 20%. On the other hand, 
experiments with additional scattering material, and comparison with the 
backscatter peak from the 661 kev. gamma ray of Cs!*7 showed that the ob- 
served peak was of the correct order of magnitude to be due to backscattering 
alone. Au’ was therefore not detected, and at worst the amount of Au!” 
present in the Au’ source was less than (0.2 K 1.3)% = 0.26%. The corre- 
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sponding half-life error is less than 0.04%. The amount of Au™ present, 
estimated from the neutron flux during irradiation, was less than 0.001%. 
The source was also searched for gamma rays higher in energy than 412 kev., 
and no radiations were found comparable with the 0.2% abundant 1.089 
Mev. gamma ray of Au" (2), 

The decay of the Au" source foil was measured using a TQQ quartz fiber 
electroscope (ionization chamber) in the manner described by Sargent, Yaffe, 
and Gray (6). The background of the instrument before the experiment 
was 0.01 div./min., and the initial activity of the sample was 3 X 10‘ times 
background. The decay was followed without removing the source from the 
electroscope until the readings returned to the previous background. Readings 
were taken once or twice daily for 10 half lives, and once daily thereafter 
until 18 half lives had elapsed (decay factor 2.5 X 10°). The readings were 
corrected for temperature and pressure variations in the usual way (6), and 
for background; no reading lasted long enough to require correction for 
duration. When plotted on semilogarithmic paper, the corrected readings 
appeared to lie exactly on a straight line over the whole time of measurement, 
except for a scatter of the last few points as the readings approached back- 
ground. 

RESULTS 


The 38 points lying in the first 10 half lives (final reading 30 times back- 
ground) were fitted to a straight line by least squares, the result being 2.6994 + 
0.0016 days for the half life. The indicated error is the standard deviation 
computed from internal consistency of the points. The upper limit of error 
in the half life due to the possible (but unlikely) presence of Au’ in the Au’ 
source is 0.04% = 0.0011 days. The result of the present measurement may 
therefore be given as 2.699 + 0.003 days, the indicated probable error being 
slightly greater than the sum of the errors just mentioned. 

While this work was in progress, Lockett and Thomas (5) reported the 
result of a measurement of this half life, using a method similar to that 
reported here but without the precautions against Au’ contamination. 
Their value is 2.697 + 0.003 days. The close agreement between the average 
of Table I, the value of Lockett and Thomas, and the value from the present 
work makes it unlikely that any ‘‘Au’® error’ exceeding 1 part in 1000 is 
present in the results. From all these results, the best value for the half life 
of Au is 2.698 + 0.003 days, the probable error being a conservative one. 


The authors wish to thank the National Research Council of Canada for a 
grant-in-aid and Atomic Energy of Canada Limited for the loan of equipment. 
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THE MOMENTUM AND VORTICITY TRANSFER THEORIES 
OF TURBULENT HEAT TRANSFER! 


By A. W. Marris? 


ABSTRACT 


The case of heat transfer by a cylindrical turbulent region is examined further 
theoretically from the standpoint of the vorticity transfer analogy. The radial 
distribution of the eddy diffusivity for vorticity is considered and, on the 
logarithmic velocity distribution law for fully developed turbulence, this quan- 
tity is found to be negative throughout an interval at the outer boundary of 
the turbulent region. When this region is excluded from the relevant integrals, 
results are obtained for the Nusselt modulus and radial temperature distribution 
for the particular case of Prandtl number equal to the ratio of the eddy diffu- 
sivities for vorticity and heat, and are compared with the corresponding results 
on the momentum transfer analogy theory. 


SUMMARY OF NOMENCLATURE 


é€y = eddy diffusivity for momentum. 
ey = eddy diffusivity for vorticity. 
€y = eddy diffusivity for heat. 
Pr = Prandtl number for fluid = cyu/s. 
c = specific heat of fluid. 
“ = molecular viscosity. 
s = molecular thermal conductivity. 
Re = Reynolds’ number = 27 oum/v. 


To To 
Um = mean flow velocity defined by um = - ur ar/ f r dr. 
0 0 


ro = outer radius of flowing fluid. 
v = kinematic molecular viscosity = u/p. 
p = density of fluid. 
Q = Rev(f/8) . 
f = Fanning friction factor. 
Nu = Nusselt modulus = 2roh/s. 
h = heat transfer coefficient of flow boundary. 
r = radial distance measured from axis of flow. 
x= r/Tr. 
u = mean flow velocity at r. 
t = fluid temperature at r. 
to = fluid temperature at 7o. 


To To 
tm = mean mixed fluid temperature defined by tm = J urt ar / fr dr. 
0 0 


rt = fluid shear at 7. 
To = fluid shear at 7o. 


u 2 : , ; 
U= u/ 4/* = #4/2 = nondimensional mean flow velocity at x. 
m 


1Manuscript received March 12, 1954. 
2Present address: Physics Department, University of British Columbia, Vancouver, B.C. 
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M, N = constants in the turbulent velocity distribution law 
U = M+ WN log Q(1 — x). 
8 = dimensionless boundary of von Karman’s buffer layer, i.e., 
x = 1 — (8/Q) signifying outer boundary of turbulent region 
for which U = M + N log Q(1 — x) holds. 
M . N log B. 
M + Niog N. 


INTRODUCTION AND SUMMARY 


On the “mixing length” theory of turbulence the rapidly moving eddies are 
assumed to transfer properties such as heat momentum or vorticity just as 
in the classical kinetic theory of gases these transferable properties are assumed 
to be transported by molecular agitation. Thus setting momentum as a 
transferable property one is lead to an additional turbulent viscosity expressed 
as an eddy diffusivity for momentum ey, and likewise for heat transfer one 
obtains the eddy diffusivity for heat eg, a turbulent addition to the molecular 
thermal diffusivity. Assuming a relationship between eg and ey one may, in 
particular cases, integrate the equations to give explicit expressions for the 
heat transfer coefficient and temperature distribution throughout the turbulent 
fluid (4). 

A necessary assumption implicit in this momentum transfer theory is that 
the fluctuating pressures that exist in the turbulent field are ineffective as far 
as the mean transport of momentum is concerned. This difficulty may be 
overcome (2, 6, 7) by taking vorticity instead as the property transferred by 
the eddies, in which case one obtains the turbulent viscosity as ey, an eddy 
diffusivity for vorticity. The assumption of a relationship between ey and ey 
yields an alternative heat transfer result on the vorticity transfer analogy. 

It has been stated that this vorticity transfer theory breaks down in the 
neighborhood of solid boundaries (3). 

In this paper it is shown that on the well-defined logarithmic turbulent 
velocity distribution law the eddy diffusivity for vorticity ey is in fact negative 
throughout a region of thickness O(1/Q) bounding the turbulent core, a 
region the width of which is of a higher order than the O(1/Q) assumed by 
von Karman to be the extent of the laminar and ‘‘buffer’” layers adjacent to 
the boundary. Here Q = Re /(f/8). By cutting this region from the range of 
integration the heat transfer to be expected by a turbulent core or wake on 
the vorticity transfer analogy is determined for the case P’ = 1, where 
P’ = (e;/ey)Pr, and compared with that on the momentum transfer analogy 
for P = 1, P = (ex/ex)Pr. 

One obtains for 1/Nuf the series 5[F,(Q)/Q*" as opposed to the series 
>G,(Q)/Q’ obtained in the momentum transfer analogy. Here F,(Q), G,(Q) 
are determined functions of Q. The term O(1/Q*) dominating, the value of 
Nuf is found to be very much less than its corresponding value on the mo- 
mentum transfer analogy and also to increase with Q much more slowly. It is 
apparent that for a given value of the Prandtl modulus the results of the two 
representations could not be made to coincide by any choice of ey/ey or 
ey/ew independent of Q. 


7 
L 
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THE EDDY DIFFUSIVITIES FOR MOMENTUM 
AND VORTICITY 
For the case of turbulent flow through a tube of circular section one has, 
on the momentum transfer theory, the hydrodynamical relationship 


(1) 1 P12 ly +2) |, 


Equation [1] expresses the mean pressure gradient in terms of the radial 
velocity gradient, the kinematic viscosity, and the eddy diffusivity for mo- 
mentum. Alternatively for a given radial velocity distribution the equation 
may be regarded as expressing €y as a function of the radius. 

Thus with the Fanning friction factor f defined by 


[2] te gl oP cs acai 
= Um dz 3pUm . 


a nondimensional velocity U by 


‘ aval areas 


a nondimensional distance from the tube axis by 
[4] x =r/to, 


equation [1] integrates to 


5] “t= — (ew +») & 


so that 
ew a 
[6] - dU /ax 


On the vorticity transfer analogy the equation corresponding to [1] would be 


1 oP 1 a ou 
[7] = (ey + nl? = (, a. 


an equation which is seen to express the eddy diffusivity for vorticity, ey, in 
terms of the second derivatives of the velocity with respect to the radius. 
One obtains, in place of [6], the result 


[8] 41m —oe/42 (<%). 

Utilizing now the experimentally well defined turbulent velocity distribution 
[9] U = M+Niog Q(1 — x) 
where M = 5.5, N = 2.5, 


equations [6] and [8] yield respectively 


[10] = = £ =i ~s) = 1, 








422 CANADIAN JOURNAL OF PHYSICS. VOL. 32 


, 2 9 
[11] 4g ae x(l1— <x)" —1. 


See Fig. 1 for the graphs of these functions. 

It is seen that while the function €,;/y is negative throughout intervals of 
length O(1/Q) at the boundary « = landatx = 0, the corresponding function 
ey/v is negative throughout an interval of length O(1/Q!) at the boundary 
x =1. 








Oo 'Q2 ‘4 (G6 O68 
x 





Fic. 1. Radial distribution of eddy diffusivities for momentum and vorticity for Q@ = 1250 
and Q = 5000. 


It is because of this latter interval at the boundary of the turbulent region 
that impossible results were obtained in the previous integrations for the 
Nusselt modulus and turbulent temperature distribution (4). The range of 
integration was the extent of the turbulent core 0 < x < 1 — (8/Q) according 
to von Karman’s buffer layer hypothesis. It is now seen that on the vorticity 
transfer theory the adjacent-to-boundary interval for which ey/v is negative 
so that ey and »v will cancel each other is of a higher order of thickness than 
8/Q, and that instead of setting 8 = 30 one would have to set 6 = RQ}. 

Values for the product Nuf and the turbulent temperature distribution on 
the vorticity transfer theory may be obtained by cutting from the range of 
integration the intervals over which ey/v is negative. The limits could be 
found exactly by solving the cubic 
[12] ze x(1—x)?-1=0, 
which has three real roots; however, because of the form of the velocity 
distribution the range of x considered as covering the turbulent region will be 
N/Q <x <1—(N/Q)?. The ratio ey/y is positive in this range, its values at 
the end points being 1 — (4N/Q) + O(1/Q") and 1 — 2(N/Q)}, values very 
small in comparison with its mean value of order Q. 
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It has already been shown (4) that to O(1/Q) the result for the temperature 
distribution and Nusselt modulus as derived on the momentum transfer anal- 
ogy is unaffected by the negative value of €y,/v in the range 0 < x < 6/Q at 
the center of the turbulent core. 

In what follows the heat transfer to a cylindrical turbulent core extending 
from x = 0 to x = 1 — (8/Q) will be considered from the standpoint of the 
momentum transfer analogy for the particular case P = (ez/ey)Pr = 1 and 
compared with that to a similar turbulent core extending from x = N/Q to 

= 1 — (N/Q)! considered from the standpoint of the vorticity transfer 
analogy for the particular case P’ = (e,/ey)Pr = 1. The same experimentally 
defined velocity distribution U = M + N log Q(1 — x) will be used in each 


case. 


EXPRESSIONS FOR THE NUSSELT MODULUS AND TEMPERATURE 
DISTRIBUTION ON THE TWO THEORIES 
It has been shown (4) that on the momentum transfer theory, for a turbulent 
core extending over the range 0 < x < 1 — (8/Q) and for a linear longitudinal 
temperature gradient, the radial temperature distribution and the product 
of the Nusselt modulus and friction factor are given by the following ex- 
pressions: 





ot _nug/f 
[13] oa Nu4 > I (x3), 
1 1—(8/Q) 
[14] Nuf = J Ux; I (x3) dx3 
where 


1—(8/Q) 
s 2 PQx2 
[15] I(x3) = J Xe {a -P)- Fon. |" rt Ux, dx dx. 
For the particular case P = (ey/ey)Pr = 1 equation [15] becomes 


1—(8/Q) Zs 
[16] I(x3) = - f a |. U x dx, dx». 


Writing y = M+NlogB 





so that U = M + N log Q(1 — x) = [y — N log (6/Q)] + N log (1 — x), 
one finds 

2 N N N 
J ‘Varden = (MH — Mog 8) 2p — May — Nc — 1%) tog (1 - =) 
so that 


sab ay “orl i(2y—N) —3N log (8/Q)}x2” 


x2 (1 — x2) 


isa iNx+}N(1— —x2") log g(l==)) 46 
x2 (1 — x2) 
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17 2y-3N WN N B 2y-N WN B\N 
17] = ——_ a log 8 log go + = a log aw log (1—xs) 
-<! log gO Hs) x3) += 20 W"| tog x3 log (1 — xs) + > ne + o(2) 
and 
[18] (x) = logx log (l — x) + SS 2 = +> a = )==>- 


Equations [13], [17], [18] give the temperature distribution once the Nusselt 
modulus is determined (4). 


For the Nusselt modulus 


1 1—(8/Q) 
Nuf = J Uxs3 I (x3) dx3 


oh res e) + N log (1 — «| x31 (x3) dxs 


2y-3N _N aie B 
-( 4 — 2) o) g'%Q 

Q24—-N N 7N — 3 3 N 
+ (25% - Zee 8) ag hae 8) 








7 — 6yN aut B 








N 
"7 x ( iO Q 
7N’y _3Ny*_ 71IN' , MS J 1 ia) 
+|4 8 329 to4elot\e 
Putting 
[20] M =5.5, N = 2.5, y = M+ Nlog30 = 14.0, 


one obtains finally, on the momentum transfer theory, 


1 47. fl 30 138 30 

Pll Ny" —~@ (ios a0) +9 (Iog1032) 
_ 15.6, 30 55.4 es ) 
“o rg 6 TON 


On the vorticity transfer analogy the equations corresponding to equations 
[13], [14], [15] will be (4) 
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[22] ns a = wag/L £ ies), 
1 1-(w/Q)4 
[23] Nu J Us J (x3) dxs 
where 
1-(N/Q)4 -1 pt 
=~ —1 — >, om 2P’Q | 7, 

[24] a J. ™ [a 3 8(x20U/ x2) / Ox2 ia See 

For the particular case P’ = (eg/ey)Pr = 1 equation [24] becomes 

1-W/@)4 4 Za 

: 8 2) 

[25] J (x3) = ~ 20 aI. Xe Ox («, Oxe/ « . dx dx» 


corresponding to equation [16]. 
Writing L = M+ Nlog N 


so that U = M+ Nlog Q(1 — x) = [L — N log (N/Q)] + N log (1 — x), 





one has 
Moe ee, ee 
ec oe — x2") log (1 — x2) +O (2) , 
so that 
ctor (HQL—W)— 4 og /Q)\xs" 
J (xs) = ob [x2(1 — x2)] 
_ 4Nx.+4N(1—x:*) log (1—x2)—O 0) ax 
[x2(1 — x2)] 
2L—3N WN N 1 
[26] = F(Q) — (22aN _ 4 q/0 2%) Gi-s 
-3 {oe—a) _ (i=ss) _ 2 | tog x3 log (l—x3s) + z fica + [log (1 ~ x)" 
where 
2L—3N WN, N\(Nn\' N ,1(/, N\?\w 
[27] F(Q) = (225-88 at log *)(%) + (10g 0 aa (tog ¥) )z 


+o(d). 


Equations [22], [26], [27], [18] give the temperature distribution once Nu 
is determined. 
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For the Nusselt modulus 


1 1—(N/Q)4 
Nuf = J Ux; J (x3) dx3 
+ /Q 





eo | (z Wie x) + N log (1 — x) | x3 J (x3) dxs 


N/Q 


(22 —~3N WN ay 
== = — = lor —) 1 — 
2 4 28 0} \o 











n 
28] -i(% son | Sto Ae res ( Nw)’ 
! “a. 6) \o) + 169 VF oO 
N?(6L — 5N) ( ay 
320 =o 
(2L —3N)(2L —N)N, N ae et 4)2 ( 1 ) 
aa 320 log +3 (i+ . + 0\ oy) - 
Putting 
[29] L=M+NlogN = 5.5 + 2.5 log 2.5 = 7.80, 


one obtains finally, on the vorticity transfer theory, 





[30] Nuf 7 0.791 (2.02 - 2.88 logw2®) 2 at = logio 4) 
35.5 ( 2. a) 19.0, 2:5 4 4.30 o( 1 ) 
a eh he Se a het = + ae 
0 O10 QO + QO “0 + 0 +0 of” 
RESULTS 


Tables I and II give the values of the Nusselt modulus as calculated from 
equations [21] and [30] and von Karman’s friction factor law 


[31] 1/f? = 4.0 logio Q + 1.41 

where 

[32] Re = 11.3 Q logis Q + 3.98 Q. 

TABLE I 
MOMENTUM TRANSFER ANALOGY CASE Pr = ey/ex 

Q Re X 10 Nuf Nu 
1000 0.379 2.19 393 
2000 0.825 2.77 588 
3000 1.30 3.30 775 
5000 2.29 4.22 1110 


20000 10.5 9.21 3180 
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TABLE II 
VoORTICITY TRANSFER ANALOGY CASE Pr = ey/ez 

Q Re X 10° Nuf Nu 
1000 0.379 0.566 102 
2000 0.825 0.634 134 
3000 1.30 0.684 160 
5000 2.29 0.756 199 
20000 10.5 1.02 354 


The variation of Nuf with Q = Re+/(f/8) according to the two representations 


is shown graphically in Fig. 2. 
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Fic. 2. Variation of Nuf with Q. 


If the values of the Nusselt modulus obtained on the momentum transfer 
analogy for a wholly turbulent core are compared with the values obtained for 
a turbulent core surrounded by a laminar flow layer for which ey = eg = 0 (4), 
it is seen that the existence of such a laminar layer reduces the Nusselt modulus 
by a factor of 2. 

It is with this in mind that one may see the reason for the low values of 
Nuf obtained on the vorticity transfer theory and for the small variation of 
Nuf with Q. From the graphs of Fig. 1 it is seen that while the maximum value 
of ey/v is always greater than the maximum of ey/v, the vorticity diffusivity 
maximum is displaced towards the tube axis and ey/v is less than ey/v through- 
out the whole region from x = 3 to the boundary. So too, eg on the vorticity 
transfer representation is less than eg on the momentum transfer theory if one 
persists in the assumption ey/e4 = €7/ey = k. Efficient turbulent heat transfer 
depends first and foremost on the relevant eddy diffusivities reaching their 
high values as near to the boundary as possible. 

If (Nu) y,(Nu) y are written for the Nusselt moduli obtainedon the momentum 
and vorticity transfer theories respectively, the variation of (Nu) s/(Nu) y with 
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Q is as shown in Fig. 3. Eckert (1) has remarked that in wakes satisfactory 
results on the momentum transfer analogy could only be obtained by putting 
€/€u = 2. In terms of the previous theory this is tantamount to saying that 
the Prandtl number must be doubled in order to predict reconcilable results, 
or that Nu is only about 3 of its theoretically predicted value (4). 

If these observations were made at Reynolds numbers of around 10,000- 
20,000, such a value for Nu would be predicted on the assumption that the 
vorticity transfer representation held for wakes with the ratio ég/ey set as 
unity (3). Further experiments for systems for which the vorticity transfer 
representation might be expected to hold, at high Reynolds numbers, would 





Q«10° 





O 2 4 6-6 WO. de le ses WG 
Fic. 3. Variation of the ratio (Nu)yw/(Nu)y with Q. 


be of extreme interest in showing whether the Nusselt modulus is in fact 
lessened from that predicted by the momentum transfer analogy by the ratio 
indicated by Fig. 3. 

Attention is drawn to certain experimental results obtained by the author 
for the heat transfer to air pressurized in order to extend the Reynolds number 
range (5). At the higher Reynolds numbers it was found that the experi- 
mentally determined values of Nuf dropped off to about one half that predicted 
by the momentum transfer analogy theory. Whereas these results are by no 
means conclusive, the experimental setup then used being subject to considera- 
ble criticism, the results obtained could be partly explained by the assumption 
that vorticity transfer dominated at the higher Reynolds numbers. It would 
be of interest to repeat this work under perhaps simpler experimental con- 
ditions. 


CONCLUSION 


From an examination of the distribution of the eddy diffusivity for vorticity 
across the radius of cylindrical turbulent ‘‘core’’ it is found that a result may 
be obtained for the turbulent heat transfer coefficient on the vorticity transfer 
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analogy theory provided that an interval of thickness O(1/Q*) at the outside 
boundary of the core, an interval through which the eddy diffusivity for 
vorticity is negative and theoretically works against the molecular viscosity, 
be excluded from the range of integration. The values of Nuf obtained for the 
particular case Pr = ey/eg on the vorticity transfer analogy are compared 
with those obtained on the momentum transfer analogy for the particular case 
Pr = €x,/ez. It is found that the former depend on Q} whereas the latter 
depend on Q. The small heat transfer coefficients obtained on the vorticity 
analogy representation result from the fact that the distance of penetration 
into the turbulent region required before ey reaches its high values is larger. 
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ON A TWO-FLUID MODEL OF MATTER! 


By F. A. KAEMPFFER 


ABSTRACT 


The theory of Pauli and Weisskopf is reformulated in hydrodynamical terms 
by introduction of a suitable density p and a suitable velocity potential ¢. The 
Hamiltonian of the system is expressed in terms of p, ¢, and their canonical 
momenta. In accordance with a general program proposed by the author the 
transition to quantum theory is carried out along the lines of quantum hydro- 
dynamics. The eigenvalues of the Hamiltonian corresponding to excitations of 
the motion with no net mass flow are obtained. Since these excitations do not 
give rise to an electromagnetic field, they are tentatively identified with neutrino- 
like particles. The picture emerging from these considerations is that of two 
interpenetrating fluids of positive and negative charge, in which different types 
of elementary particles appear as different types of excitations of the motion. 


The most remarkable feature of all known forms of matter is the experi- 
mentally established fact that any electric charge density carried by matter 
can be either positive or negative, while any mass density carried by the 
same matter is always positive only. 

A well-known example of a description that recognizes this fact is the 
theory of Pauli and Weisskopf (5), in which a complex scalar field y, y*, 
satisfying the classical Klein-Gordon equations 


(1) (O*—- x) = 0, (0? — x7)" = 0, 


is used to describe charged matter, with the aim of obtaining upon quantization 
of this theory a description in terms of particles of positive mass carrying a 
charge that may be either positive or negative. The quantity A is a character- 
istic length of the field, and is interpreted as the Compton wave length of the 
particles that appear upon quantization of the theory. 


The quantities 
* 
«ih <a) 
i) Pee (y dx, v Ox, 


are identified with the components of the four-vector of electric charge density, 
because they can be either positive or negative, and because as a consequence 
of the field equations [1] they satisfy the law of conservation 


[3] > <0. 


7 
The factor ¢ is also a characteristic parameter of the field and is given the 
dimension of a charge density. 
Since the equations [1] are derivable from a variational principle character- 
ized by the Lagrangian density (7) 


* 
[4] L=-¢ (x ~ “ + yy) ; 


v 


1 Manuscript received February 18, 1954. , 
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one can define as usual a stress tensor 


* * 
[5] z. a ce (x oy $ oy ay) + Liu» 


Ox, OX, CXp Ch 











such that the mechanical behavior of the substance under consideration is 
completely determined by the conservation laws 

Cy Po 
[6] x = ~ ° 
with the understanding that G,=(1/ic)T«,., Pik = Tix, Se=(C/t) Tis, H= —T 4 
are interpreted respectively as momentum density, mechanical stress density, 
energy flux density, energy density of the field. 

It was the purpose of the work reported in this paper to show that the theory 
of Pauli and Weisskopf can be interpreted in terms of a two-fluid model, in 
order to obtain a formulation in hydrodynamical terms which, according to 
a general proposal by the author (4), should then be subjected to a quantization 
procedure along the lines of quantum hydrodynamics in the hope that a 
unified theory of elementary particles in which different types of particles 
appear as different quantized excitations of the motion of the field under 
consideration might thus be obtained. 

Considering the wave functions y, ¥* as dimensionless quantities we shall 
now describe the Pauli-Weisskopf theory in terms of a mass density 


[7] p= wy 

and a velocity four-vector 

[8] u,=v,+ 4, 

where 

[9a] v, = rc = (In ¥), 
Ox, 

[ 96] w, = de ~ (in y*). 


The factor » is assumed to be positive and is given the dimension of a mass 
density. We are thus assured that the mass density associated with the field 
will always be positive. The condition that the electric current density four- 
vector should be capable of describing positive and negative charge associated 
with the field and should fulfill the law of conservation of charge can be satis- 
fied by defining this four-vector similar to [2] by 
* 
[10] Ss, = — tedc (v 6 2) ; 


Ox, 


Expressed in terms of p, v,, w, this four-vector takes the form 


1€ 
11 = os v)- 
[11] s Oe Wy) 


Considering € as a positive factor we see thus that the velocity v, describes the 
motion of negative charge, while w, describes the motion of positive charge. 
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In view of the definitions [9] we may therefore also say that In y is the velocity 
potential determining the flow of negative charge, while In y* is the velocity 
potential determining the flow of positive charge. It should also be noted that 
a net flow u, can happen only if there is a gradient in the mass density present:T 


Ac A 
[12] vy + w, = de 5 fin (y"y)} = oF 
Conversely, a flow v, = — w,, which gives rise to an electric current, is not 


accompanied by any net flow of mass. Equation [12] plays the part of an 
auxiliary condition imposed upon the velocities. It can be used to eliminate 
either v, or wy. 

Introducing the velocity potential ¢ by 


pets ie de a 
[13] v, — Wy = dE ee {in (+ ae 
we have 
1 de (op , oe) 
{140} es . (2243 ax,/]’ 
_1d¢ 20 
[14d] Ww =>5 : (# 1 ax, 


and can thus describe the Pauli-Weisskopf field in terms of the two field 
functions p and ¢ and their derivatives. The Lagrangian density [4] takes now 
the form (after multiplication with \?4 for dimensional reasons) 





vee ae . dp ap , a6 at.) 
[15] ee 2 edo (2 ax, * Ox, dx.) ° 
This permits introduction of canonical momenta defined as usual by 
1 aL eo 
[16a] "eo ic A(Op/dx,)  2ue °” 
1 aL - 
[165] ™ = 5 3(a¢/dm)  Qnt 
so that 
2 
[17a] p= ss, 
A" p 
-_ . Qu? 1 
[17] o= ru . Ty. 


We obtain thus a Hamiltonian density H = 2,9 + mo — L of the Pauli- 
Weisskopf field in terms of the density p and velocity potential ¢ and their 
canonical momenta: 


tWe imagine a In p expanded in — = p—u. If one does not wish to make this 


approximation, don po de, has to be replaced i? on ~ in subsequent equations. 
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2 
[18] H= prt fe (x, + ++ Qs r pf(grad p)* + (grad ¢)'}. 


According to the program proposed by the author (4) the transition to 
quantum theory should now be carried out by subjecting the field variables 
and their momenta to commutation relations 


[19] [o(x); mo(x’)] = thd(x — x’), [o(x); wp(x’)] = thd(x — x’), 


with all other commutators vanishing, in analogy to the commutation relations 
of quantum hydrodynamics (8). 

Complete diagonalization of the resulting Hamiltonian is certainly very 
difficult, and will not be attempted here. However, one can easily obtain the 
eigenvalues of the Hamiltonian for those excitations of the motion that are 
not accompanied by any net flow of mass, i.e. for p = po = constant. In this 
case we have, because of equation [17a], 7, = 0, and the Hamiltonian density 
reduces to 


kee 


[20] H® = poc® +752 - T + (grad ¢)”. 


Upon resolution of the operators into their Fourier components 


[21] ¢=VI> Khe, m= re P,e™", 
k 
we obtain as total Hamiltonian of the system 
[22] 9° = fa°av = po?’V + » (4 PyP_.+ —S a oh FP). 


This may be further simplified by introduction of the creation and annihilation 
operators . a, 


Ad“ cpok 
[23] .- / 36 (auton), Pr= ig/ Bae (ak — a_x) 
giving 
(24) ©° = pc?V + Hic) kay an + axax) 


with eigenvalues 
[25] E® = po? V + hic > k(Nut 4) (Ny positive integers). 
k 


Following a suggestion by Dirac (2), which was later supplemented by 
the author (4), we should now interpret 


«we 
[26] fae '3 &* eu ” ax, 





as the electromagnetic potentials. This identification is consistent with the 
present theory, because equation [3] expresses then the gauge invariance of 
the potentials. Consequently, the electromagnetic field tensor produced by 
any excitation should take the form 
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2dr _ dn _ 2e {ae a¢ ap a) 


[27] oo = dX, ax, = eu \Ox, OX, OX, Ox, 


Ox, OX, OX, Ox,/° 


It is apparent, with this interpretation, that excitations corresponding to 
p = po = constant do not give rise to an electromagnetic field. We conclude 
therefore, tentatively, that the excitations of the motion which contribute to 
the energy of the system the eigenvalues [25] are to be identified with neutrino- 
like particles. 

Success or failure of this approach rests now on the outcome of the complete 
diagonalization of the Hamiltonian resulting from [18]. It is clear that exci- 
tations corresponding to variable p may give rise to electromagnetic fields 
according to equation [27], and may therefore be identified with photonlike 
and electronlike particles, depending on whether the expectation value of 
div E in these states of the Pauli-Weisskopf field vanishes or not. 

It might also be worth while to look into the possibility of generalizing the 
current density four-vector s, so that it will include vortex flow of the velocity 
fields by the introduction of two Clebsch potentials é, 7 


__ ae, (ag =) 
aa ee oe (3 tea, 


as has been suggested by Dirac (3) and Takabayasi (6). It is conceivable that 
this might lead to a description of the spin of the excitons of the Pauli- 
Weisskopf field in terms of and y. This has already been suggested by 
Bohm (1). 

The author submits these considerations as a possible starting point for a 
unified theory of elementary particles. 
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